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Abstract 

One of the most efficient methods to obtain the vacuum expectation values for the physical 
observables in the Casimir effect is based on using the Abel-Plana summation formula. This 
allows to derive the regularized quantities by manifestly cutoff independent way and to present 
them in the form of strongly convergent integrals. However the application of Abel-Plana formula 
in usual form is restricted by simple geometries when the eigenmodes have a simple dependence 
on quantum numbers. The author generalized the Abel-Plana formula which essentially enlarges 
its application range. Based on this generalization, formulae have been obtained for various types 
of series over the zeros of some combinations of Bessel functions and for integrals involving these 
functions. It have been shown that these results generalize the special cases existing in literature. 
Further the derived summation formulae have been used to summarize series arising in the 
mode summation approach to the Casimir effect for spherically and cylindrically symmetric 
boundaries. This allows to extract the divergent parts from the vacuum expectation values for 
the local physical observables in the manifestly cutoff independent way. Present paper reviews 
these results. Some new considerations are added as well. 
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1 Introduction 



The Casimir effect is among the most interesting consequences of quantum field theory and is es- 
sentially the only macroscopic manifestation of the nontrivial properties of the physical vacuum. 
These properties may be determined from the responce of the vacuum state to classical external 
fields or constraints. The simplest case is realized by boundary conditions on quantized fields. 
Such conditions modify the zero point mode spectrum and as a result can change the energy 
of the vacuum. This change is manifest as an observable Casimir energy. Since the original 
work by Casimir in 1948 |jj many theoretical and experimental works have been done on this 
problem, including various types of boundary geometry and non-zero temperature effects (see, 
e.g., [||, ||, |], H and references therein). Many different approaches have been used: mode sum- 
mation method, Green function formalism, multiple scattering expansions, heat-kernel series, 
zeta function regularization technique, etc.. From the general theoretical point of view the main 
point here is the unique separation and subsequent removing of the divergences. Within the 
framework of the mode summation method in calculations of the expectation values for physical 
observables, such as energy-momentum tensor, one often needs to sum over the values of a certain 
function at integer points, and then subtract the corresponding quantity for unbounded space 
(usually presented in terms of integrals). Practically, the sum and integral, taken separately, 
diverge and some physically motivated procedure, to handle finite result, is needed. For a num- 
ber of geometries one of the most convenient methods to obtain such regularized values of the 
mode sums is based on the using of the Abel-Plana formula (APF) §, g, §. In @ this formula 
have been used to regularize scalar field energy momentum-tensor on backgrounds of various 
Friedmann cosmological models. Further applications to the Casimir effect for flat boundary 
geometries with corresponding references can be found in Abel-Plana formula allows (i) 
to extract by cutoff independent way the Minkowski vacuum part and (ii) to obtain for the 
regularized part strongly convergent integrals, useful, in particular, for numerical calculations. 
However the applications of APF in usual form is restricted by the flat boundary cases when 
the eigenmodes have simple dependence on quantum numbers. 

In [10, 11] the APF was generalized (see also JT^j). The generalized version contains two 
meromorphic functions. Choosing one of these functions in specific form APF in usual form is 
obtained. By applying the generalized formula to Bessel functions in [jll], 12] summation formulae 
are obtained over the zeros of various combinations of these functions. In particular, formulae 
for Fourier-Bessel and Dini series are derived. From these formulae by specifying the constants 
and choosing the order of Bessel function equal to 1/2 one obtains a simple generalization of 
APF for the case of a function having poles. It have been shown that from generalized formula 
interseting results can be derived for infinite integrals involving Bessel functions. Further the 
obtained summation formulae are applied to regularize the vacuum expectation values for the 
energy-momentum tensor components of the electromagnetic field in the Casimir effect with 



spherically [|1J3, 14, 15, 1£] and cylindrically symmetric [17, 18] boundaries. As in the case of 



flat boundaries the using of generalised Abel-Plana formula allow to extract in manifestly cutoff 
independent way the contribution of the unbounded space and to present regularized values in 
terms of exponentially converging integrals. 

The present paper reviews these results and is organized as follows. In section 2 the general- 
ized Abel-Plana formula is derived and as a special case usual APF is obtained. It is indicated 
how to generalize this formula for the functions having poles. The applications of generalized 
formula to Bessel functions are considered in the next section. We derive two formulae for the 
sums over zeros of AJ v {z) + BzJ' v {z). Specific examples of applications of the general formulae 
are considered. For v = 1/2, B = and for analytic function f(z) the APF is obtained. In 
section 4 from generalized Abel-Plana formula by special choice of function g(z) summation 
formulae are derived for the series over zeros of the function J u (z)Y u (Xz) — J u (\z)Y u (z) and 
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similar combinations with Bessel functions derivatives. Special examples are considered. The 
applications to the integrals involving Bessel functions and some their combinations are dis- 
cussed in section 5. A number of interesting results for these integrals are presented. Specific 
examples of applying these general formulae are described in the next section. In section 7 by 
using generalized Abel-Plana formula two theorems are proved for the integrals involving the 
function J l/ (z)Y fJ/ (Xz) — J fl (Xz)Y I/ (z) and their applications are considered. The following sec- 
tions are devoted to the applications of generalized formula for the calculations of the regularized 
vacuum expectation values of the electromagnetic energy-momentum tensor inside (section 8) 
and outside (section 9) a perfectly conducting spherical shell, and for the region between two 
perfectly conducting spherical surfaces (section 10). In sections 11-13 the similar problems for 
the cylindrical surfaces are considered. The section 14 concludes the main results considered in 
this paper. 



2 Generalized Abel-Plana formula 



Let f(z) and g(z) be meromorphic functions for a < x < b in the complex plane z = x + iy. Let 
us note by Zf^ and z g ^ the poles of f(z) and g(z) in region a < x < b, respectively. Assume 
that Imzf^ 7^ (see however the Remark to Lemma). 



Lemma. If functions f (z) and g(z) satisfy condition 

rbznih 

lim / [g(z)±f(z)]dz = 0, 
then the following formula takes place 

r-b 



[° f{x)dx = R[f(z),g(z)} - \ P°° [g(u)+sgn(lmz)f\u)r u Z b a X z z dz, (2.2) 

Ja ^ J —ioo 



where 



R[f(z),g(z)]=7ri 



ReS z=z 9 , k g( Z ) + Res lmz Lk >of(z) ~ Res Im 2/ , fc <0/(^) 
. k k k 



(2.1) 



(2.3) 



Proof. Let us consider a rectangle Ch with vertices a±ih, b±ih described in the positive sense. 
In accordance to the residue theorem 



/ g(z)dz = 2-ni ^ Res z=Zg k g{z), 



(2.4) 



where rhs contains sum over poles within Ch- Let and C h denote the upper and lower halfs 
of this contour. Then one has 

/ g{z)dz = [ [g(z) + f(z)]dz + / _ [g(z) - f{z)]dz - [ f{z)dz + / _ f(z)dz. (2.5) 
JC h JC h JC h JC h JC h 

By the same residue theorem 

/ f(z)dz - ( f(z)dz = 2 I" f(x)dx + 2vri 



J2Resi mZf k<0 f(z) - ^2Resi mZf k>0 f(z) 
k k 



(2.6) 
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Then 



[g(z) ± f{z)}dz 



±ih 



[g(u)±f(u)r u Z b a + + z z dzT 



b±ih 



azrih 



[g(z) ± f{z)]dz. (2.7) 



Combining these results and allowing in ( |2.4| ) h —* oo one obtains the formula Q2.2| ). ■ 

If the functions f(z) and g{z) have poles with Recife = a,b (i = f,g) the contour have to 
pass round these points on the right or left, correspondingly. 

Remark. The formula ( |2.2j ) is valid also when the function f(z) has real poles Zf n , Imzj. ^ = 
in the region a < Hez < b if the main part of its Laurent expansion near of these poles does 
not contain even powers of z — z^ n . In this case on the left of the formula ( |2.2[ ) the integral is 
meant in the sense of the principal value, which exists as a consequence of the abovementioned 
condition. For brevity let us consider the case of a single pole z = Zq. One has 



f(z)dz 



f(z)dz 



ZQ-P 



f(z)dz + 



20 +p 



f(z)dz 



+ 



+ 2tti 



Resim^ th <of(z) - E Res i 



mz/,fc>0 



f(z) 



+ f + f(z)dz+ f _f(z)dz, (2i 



with contours and T p being the upper and lower circular arcs (with center at z = zq) joining 
the points zq — p and zq + p. By taking into account that for odd negative I 



r+ 
1 p 



(z - z ) l dz + l (z - z ) l dz 
in the limit p — > we obtain the required result. ■ 



(2.9) 



In the following on the left of (|2.2| ) we will write p. v. J a f(x)dx, assuming that this integral 
converges in the sense of the principal value. As a direct consequence of Lemma one obtains 

a 

Theorem 1. If in addition to the conditions of Lemma one has 

rb±ioD 



lim 

b—*oo Jfy 



[g{z] ± /(*)] dz = 0, 



then 



b— >oo 



lim^p.v. / f(x)dx-R[f(z),g(z)] 



a+ioo 



[g(z) + sgn(Imz)/(z)] dz, 



(2.10) 



(2.11) 



where on the left R[f(z),g(z)] is defined as (2.5), a < Rez/^, Rez 9) fc < b, and summation goes 
over poles ztj. and z g ^ arranged in order Rez^fc < Re^fc+i, i = f,g. 



Proof. To proof it is sufficient to insert in the general formula (|2.2j ) b — > oo and to use the 
condition ( 2.10[ ). The order of summation in R[f(z),g(z)] is determined by the choice of the 
integration contour Ch and by limiting transition b — > oo. ■ 

We will call the formula ( 2.1l| ) as Generalized Abel-Plana Formula (GAPF) as for 
b = n + a, < a < 1, g(z) = —if(z) cot -kz and analytic functions f(z) from ( |2.11| ) follows the 
Abel-Plana formula (APF) [|, 0, § 



lim 

n— >oo 



n+a 



f(x)dx 



1 

27 



f(z) (cot nz — i)dz 



1 

27 



a+ioo 



f(z) (cot ttz + i)dz. 

(2.12) 
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The useful form of ( |2.12|) may be obtained performing the limit a — * 0. By taking into account 
that the point z = is a pole for integrands and therefore have to be around by arcs of the 
small circle C p on the right and performing p — > one obtains 



£/(") 

n=0 



f( x )dx+±f(0) + i 



f{ix) - f(-ix) 



o2lTX Y 



dx. 



Note that now the condition (2.1) is satisfied if 

lim e- 2 ^ y \\f{x + iy)\ = 



(2.13) 



(2.14) 



uniformly in any finite interval of x. The ( [2.13 ) is the most frequently used form of APF in its 
physical applications. Another useful form (in particular for fermionic field calculations) to sum 
over the values of an analytic function at half of an odd integer points can be obtained from 
©li 

f(ix) - f(-ix) 



£/(n+l/2) = 

n=0 

By adding to the rhs of (2.13) the term 



f{x)dx 



2nx 



+ 1 



-dx 



(2.15) 



IT I 



I £ Resi m2/ tk> of(z) - £ Resj mZf k<0 f(z) - i £ Res z=Zf k [f(z) cot irz] \ (2.16) 

I k k k J 



the APF may be generalized for the case when the function f(z) has poles Zfk, ~Rezf t k > 0, 
Zf,k / 1,2,.... 

As a next consequence of ( 2.11 ) a summation formula can be obtained over the points 
z n , Rez n > at which the analytic function s(z) takes integer values, s{z n ) is an integer, and 
s'(z n ) 7^ 0. Taking in ( |2.11 ) g(z) = —if(z) cot7rs(z) one obtains the following formula [ 10 1 



s'(z n ) 



w + 



f(x)dx + 



f(ix) 



f{-ix) 



g—2iris(ix) j ^2Tris(—ix) ^ 



dx, 



where 



w 



0, if 8(0) + 0,±1,±2,... 

/(0)/[2s'(0)], if a(0) = 0,±l,±2,... 



(2.17) 



(2.18) 



For s(z) = z we return to APF in usual form. An example of applications of this formula to the 
Casimir effect is given in [10]. 



3 Applications to Bessel functions 

The formula (|2.11| ) contains two meromorphic functions and is too general. To obtain more 
special consequences we have to specify the one of them. As we have seen in previous section the 
one of the possible ways leads to APF. Here we will consider another choices of the function g(z) 
and will obtain useful formulae for the sums over zeros of Bessel function and their combinations, 
as well as some formulae for integrals involving these functions. 

First of all to simplify the formulae let us introduce the notation 

F(z) = AF(z) + BzF'(z) (3.1) 

for a given function F(z), where the prime denotes derivative with respect to the argument of 
function, A and B are constants. As a function g(z) in GAPF let us choose 
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where J u (z) and Y v (z) are Bessel functions of the first and second (Neumann function) kind. 
For the sum and difference on the right of ( |2.11| ) one obtains 

f(z)-(-l) k g(z) = -^-f(z), k = 1,2 (3.3) 



with h1 1] and h1 2) being Bessel functions of the third kind or Hankel functions. For such a 
choice the integrals ( |2.1[) and ( p.lOj) can be estimated by using the asymptotic formulae for 
Bessel functions for fixed v and \z\ — > oo (see, for example, pOj, It can be easily seen 



that conditions fl2.l|) and ( 2.10 ) are satisfied if the function f(z) is restricted by the one of the 
following constraints 

Me 2 ^ 



\f(z)\<e(x)e c ^ 

where c < 2, a > 1 and e(x) 
Bessel functions it follows that 



or 



x + ly, 



oo, 



(3.4) 



I 



b±ih 

aznih 
bznioo 



[g(z)±f(z)}dz 
[g(z)±f(z)]dz 



for x — > oo. Indeed, from the asymptotic expressions for 

* Rj>-' 2 \x±ih) 





/ 




J a 




I 

Jo 



J u (x ± ih) 



-f(x ± ih)dx 



H ( v' 2 \b±ix) 
J v {b ± zx) 



/(& ± ix)dx 



< 



< 



M ie (c -^ h 
M[/h a 

N l£ (b) 
N[/b a - 1 



(3.5) 
(3.6) 



with constants Mi, M{, iVi, iV{, and i?£ (H^) corresponds to the upper (lower) sign. 

Let us denote by \ Vjk ^ 0, k = 1, 2, 3 . . . the zeros of Jj,(z) in the right half-plane, arranged in 
ascending order of the real part, ReA^fe < ReA^fc+i, (if some of these zeros lie on the imaginary 
axis we will take only zeros with positive imaginary part). All these zeros are simple. Note that 
for real v > — 1 the function J v (z) has only real zeros, except the case A/B + v < when there 
are two purely imaginary zeros ]/], [2C[| . By using the Wronskian W[J 1/ (z), Y v (z)] = 2/ttz for (2.3) 
one finds 

R[f(z),g(z)} = 2 J2 T v {X v , k )f{X v , k ) + r lu [f(z% 

k 

where we have introduced the notations 

z 



Tu(z) 
riuim] 



(z 2 -v 2 )JZ(z)+z 2 J' v 2 (z) 
iti Resi mZfc>0 /(^) j , 

k J v\Z) 



(3.7) 
(3.8) 



TTl 



^Resi 



mzk<0 



/(*) 



0\z) 

Ju(z) 



-ir^2Resi mZk=0 f(z 



Ju{z) 



(3.9) 



Here z k Xu,i) are the poles for the function f(z) in the region Rez > a > 0. Substituting 
( |3.7D into ( 2.11] ) we obtain that for the function f(z) meroporphic in the half-plane Rez > a and 
satisfying the condition (|3.4[) the following formula takes place 



lim 2 T u {\ u , k )f(\ U)k ) + r lu [f{z)} - p.v. / 
~ +oc { k=m Ja 



a+ica M 1] {z) , i 

J^(z) 2 



/(*' 



M 2) (z) 



f(x)dx 



dz. 



(3.10) 



where ReA Ujm _i < a < ReA^™, ReA^n < b < ReA^n+i, a < Rez^ < b. We will apply this 
formula to the function f(z) meromorphic in the half-plane Rez > taking a — > 0. Let us 
consider separately two cases. 



6 



3.1 Case (a) 

Let /(z) have no poles on the imaginary axis, except possibly at z = 0, and 

f(ze™) = -e 2 ™f(z)+o(zP»), z^O (3.11) 

(this condition is trivially satisfied for the function /(z) = o(z^)), with 

_ | 2|Re^| - 1 for integer v _ _ 

1 Rev + |Re^| — 1 for noninteger v 

Under this condition for values v for which J v {z) have no purely imaginary zeros the rhs of 
Eq,(|3.10|) in the limit a — > can be presented in the form 



i r°°K v ( x ) r„-^ /(xe w2 ) + e ^ /(xe -w2)| d x+ f f ( z) i^M dz _ f f ( z) i^M dZ: 



itjp I v (x) L" • " 'i 'Jj+ J ^ / J u (z) ' J j- J u {z) 

(3-13) 

with 7+ and 7 p being upper and lower halfs of the semicircle in the right half-plane with radius 
p and with center at point z = 0, described in the positive sense with respect to this point. In 
( p,13| ) we have introduced modified Bessel functions I v (z) and K u {z) |21[ |. It follows from (3.11) 
that for z — > 

^i^f(z) = H f ( ze ~ m ) f{ze - n + o(z -i). (3.14) 
J v {z) J v {ze- m ) Jy 1 v ' y 1 

From here for p — > one finds 

D v = \ f(z) - dz- f(z) - dz = -7rRes 2=0 /(z) T — (3.15) 

Indeed, 

/• M^iz). f _ wi M 2 \ze-™), f s hP(z), , 

-D^ = / /(z) - dz -\- / /(ze ") - / r-^cfe = / /(z) - , \ dz + 



+ 



o(z- 1 )dz = if f(z)^ldz+ [ oiz-^dz, (3.16) 



7i " " Ap + +7i ^(*) 



where 7^ (7jl> see below) is the upper (lower) half of the semicircle with radius p in the left 
half-plane with center at z = (described in the positive sense). In the last equality we have 
used the condition that integral p. v. Jq f(x)dx converges at lower limit. By similar way it can 
be seen that 

D v = i[ f(z)?Mdz+ [ o(z- l )dz. (3.17) 

J ~ip +7i p Ju\ z ) ->l lp 

Combining the last two results we obtain ( 3.15[ ) in the limit p — > 0. By using (glOb , fljPD and 
( gig ) we have Q: 

Theorem 2. If f(z) is a single valued analytic function in the half-plane Rez > (with possible 
branch point at z = 0) except the poles z& (^ \ u ,i), Rezfc > (for the case of function f(z) 
having purely imaginary poles see Remark after Theorem 3), and satisfy conditions (\3.4\) and 



( 3.11 ), then in the case of v for which the function J v {z) has no purely imaginary zeros, the 
following formula is valid 



b 



lim I 2 T T v {\ V)k )f(\ U)k ) + r lu [f{z)) - p.v. t f(x)dx 1 
^+°° I k=i Jo ) 



^Res z=0 /(z)^4 " ~ T7T l^fi™*" 2 ) + e^f(xe~ m H dx, (3.18) 
2 J v (z) IT Jo l v {x) L J 
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where on the left ReA„ in < b < Re\ vn +i, < Rezk < b, and T u {\ uk ) and r iu[f(z)] are 
determined by relations (3.t) and (3.9). 

Under the condition ( |3.11| ) the integral on the right converges at lower limit. Recall that we 
assume the existence of the integral on the left as well (see section 2). The formula ( 3.18| ) and 
analog ones given below are especially useful for numerical calculations of the sums over as 
under the first conditions in ( |3.4| ) the integral on the right converges exponentially fast at the 
upper limit. 

Remark. Deriving the formula ( |3,18 ) we have assumed that the function f(z) is meromorphic 
in the half-plane Rez > (except possibly at z = 0). However this formula is valid also for some 
functions having branch points on the imaginary axis, for example, 



f(z) = h(z)l[(z 2 + ci 



±1/2 



(3.19) 



1=1 



with meromorphic function fi(z). The proof for ( |3.18| ) in this case is similar to the given above 
with difference that branch points rticj have to be around on the right along contours with small 
radii. In view of the further applications to the Casimir effect (see below) let us consider the 
case k = 1. By taking into account that 



1/2 





+ c 2 | 


Is 


+ c 2 | 




+ c 2 | 



if \z\ < c 

,2 1 1/2 e tn/2 jf j mz > c 



(3.20) 



from ( ^18|) one obtains 

lim 2 T„(\ v ,k)f(K,k) + nu[f(z)} - p.v. / f(x)da 
[e-"7i(xe n/2 ) + e^hixe-™/ 2 )] (c 2 - 
""ViCxe**/ 2 ) - e u ^h{xe-^ 2 )] (x 2 - 



1 



^TJ ft \ Y ^ Z ) 

2 JJz) 



IT Jo I u (x) 

K v (x) 



±1/2 



dx =F 



Iuix) 



±1/2 



dx, 



(3.21) 



where f(z) = f\{z) (z 2 + c 2 ) ±1//2 , c > 0. In Section 11 we apply this formula with analytic 
function fi(z) to derive the expressions for the regularized values of the energy-momentum 
tensor components in the region inside the perfectly conducting cylindrical shell. ■ 



For an analitic function f(z) the formula ( p.l8| ) yields 



y 

ti(\l k -^ 2 ) J 2 {K, k ) + \l k J' 2 {K,h) 



roc n 

/ f(x)dx + -Res z=0 .f(£ 
JO 2 



Y„(z) 
Jv{z) 



1 

vr Jo 



K V {X) 

T»{x) 



f{xe ni/2 ) + e™f(xe-™'' z ) 



7^^/2^ 



dx. 



(3.22) 



By taking in this formula v = 1/2, A = 1, B = (see the notation (|3.1| )) as a particular case we 
immediately receive the APF in the form ( ^Wj ). In like manner substituting v = 1/2, A = 1, 
B = 2 we obtain APF in the form ( |2.15| ). Consequently the formula ( 3.18 ) is a generalization 
of APF for general v (with restrictions given above) and for functions f(z) having poles in the 
right half-plane. 
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Having in mind the further applications to the Casimir effect in Sections 8 and 11 let us 
choose in ( |3.22| ) 

f{z) = F(z)J 2 +m (zt), t>0, Rev > (3.23) 
with m being an integer. Now the conditions fl3.4j ) formulated in terms of F(z) are in form 



(3.24) 



jur 2(l-t)\y\ 

\F{z)\ < \z\ee {c - 2tM or \f(z)\ < , z = x + iy, \z\ oo 



with the same notations as in fl3.4|). In like manner from the condition ( |3. 11 ) for F(z) one has 



F(ze T ' 



-F(z) + o(z 



-2m-l\ 



z -> 0. 



(3.25) 



Now as a consequence of ( p. 22 ) we obtain that if the conditions ( |3.24| ) and ( |3.25| ) are satisfied, 
then for the function F{z) analytic in the right half-plane, the following formula takes place 



00 roo 



1 



K v (x) 



- / t t \ v+r, 

for Ref > and Hev + m > 0. 



I 2 u+m (xt) \F{xe m l 2 ) + F(xe 



v+m 



{xt)dx 



dx 



(3.26) 



3.2 Case (b) 

Let f(z) be a function satisfying the condition 



f(xe 



2um 



f(xe 



(3.27) 



for real x. It is clear that if f{z) have purely imaginary poles, then they are complex conjugate: 
^Wki Vk > 0. By ( [3.27 ) the rhs of ( p. 10 ) for a — > and argA^jt = 7r/2 may be written as 



+ 



V 



E 



C P {a k ) ) Ju(z) 



f(z)dz 



+ 



E 



<?k=—iyk!—^u,k 



0\z) 

Ju(z) 



f(z)dz, 



(3.28) 

where C p {ak) denotes the right half of the circle with center at the point a k and radius p, 
described in the positive sense, and the contours 7^ are the same as in (|3.13j ). We have used the 
fact the purely imaginary zeros of J v (z) are complex conjugate numbers, as J u (ze m ) = e vm J v (z). 
We have used also the fact that on the right of ( 3.10|) the integrals (with a = 0) along straight 
segments of the upper and lower imaginary semiaxes are canceled, as in accordance of ( 3.27] ) for 
argz = 7r/2 



(3.29) 



J v {z) J v ' J u {ze-™) 

Let us show that from ( |3,29| ) for zq = x§e m l 2 it follows that this relation is valid for any z in a 
small enough region including this point. Namely, as the function f (z) Hl P \z) / J u (z) , p = 1, 2 is 
meromorphic near the point (— l) p+1 xoe m ^ 2 , there exists a neighbourhood of this point where 
this function is presented Laurent expansion 



0\z) 
JJz) 



/(*)= E 



n=—no 



[z- (-l)P+ix e^/ 2 ] n ' 



(3.30) 
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From (|3~29|) for z = xe m l 2 one concludes 

' ' i t < 

n=— no 



(X - X ) 



E 

n=— no 



(-lj^^e-™/ 2 

(X - X ) 



(3.31) 



and hence an = (— l) n On^- Our statement follows directly from here. By this it can be seen 
that 



0\z) 
CJa k ) Jv(z) 



f(z)dz 



0\z) 

C p (-a k ) Ju{z) 



f(z)dz = 2-7riRes 



z=a k 



Ju(z) 



(3.32) 



where o~k = iyk, Xu,k, argA^fc = tt/2. Now by taking into account (|3.15|) and letting p — » we 
get 10,0: ' 

Theorem 3. Let f(z) be meromorphic function in the half-plane Rez > (except possibly at 
z = 0) with poles Zk, Rez^ > and ±iyk, yk > 0, k = 1,2,... (jt X u ,p)- If this function satisfy 
the conditions (3^) and ( 3.21 ) then 

lim 2V T v (\„, k )f(\ v , k ) + r lu [f{z)\ - p.v. / 
[ £ Jo 



f(x)dx 



T E (2-5o„ fe )Res 2= „J(z) fl; ' ,!:) 



J viz) 



where on the left < Rez^ < b, ReA^n < b < ReA^+i and r\ v is defined by (3.&). 



(3.33) 



Note that the residue terms in ( |3.32f) with = A^fc, argA^. = tt/2 are equal to 4T I/ (A I/i / c )/(A I /,fc) 
and are included in the first sum on the left of (3.33). 

Remark. Let ±iyk, yk > and ±\ u k, argAj,^ = tt/2 are purely imaginary poles of function 
f(z) and purely imaginary zeros of J u (z), correspondingly. Let function f(z) satisfy condition 



/(*) 



2um 



Kze 



(( 



+ l(Z - CTfcl 



Cfc, 



0~k = i]Jk, K,k- 



(3.34) 



Now in the limit a — > the rhs of (3.10) can be presented in the form ( 3.28| ) plus integrals along 
the straight segments of the imaginary axis between the poles. Using the arguments similar 
those given above we obtain the relation ( |3.32| ) with additional contribution from the last term 
on the right of (3.34) in the form f c t- ah ) (( z — o~k)~ 1 )dz. In the limit p — ► the latter vanishes 
and sum of the integrals along the straight segments of the imaginary axis gives the principal 
value of the integral on the right of ( |3.18| ). As a result the formula ( |3.18|) can be generalized 
for functions having purely imaginary poles and satisfying condition fl3.34|) writing instead of 
residue term on the right the sum of residues from the right of ( |3.33| ) and taking the principle 
value of the integral on the right. The latter exists due to the condition ( 3.34| ). ■ 

It follows from ( |3.33 ) an interesting result. Let A^ be zeros of the function AiJ^(z) + 
BizJ'^{z) with some real constants A\ and B\. Let f (z) be an analytic function in the right 
half-plane satisfying condition (3.27) and f(z) = o(z@) for z — > 0, where (3 = max(/3„, f3 u ) (the 
definition (3 V see ( |3.12|) ). For this function from (|3.33| ) we get 



E T M»l)f(^l) = E Tu(Xu,k)f(K,k), p = v + m. 

k=l k=l 



(3.35) 
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For the case of Fourier-Bessel and Dini series this result is given in |2(J . 

Let us consider some applications of the formula ( p. 33 ) to the special types of series. Firstly 
we choose in this formula 

f(z) = F 1 (z)J^(zt), t>0, (3.36) 
where the function F 2 (z) is meromorphic on the right half-plane and satisfy conditions 

\Fi(z)\ <ei(x)e (c - t)l3/| or ^(z)] < M\z\- ai e (2 ~ t)M , \z\ -»• oo, (3.37) 

with c < 2, a.\ > 1/2, e%(x) = o(y/x) for x — > +oo, and condition 

Fi(xe™ /2 ) = -e^-^F^xe-™/ 2 ). (3.38) 

From ( 3.37] ) it follows that the integral p.v. Jq°° F\(x) J^{xt)dx converges at the upper limit and 
hence in this case the formula ( |3,33| ) may be written in the form 

°° 1 f°° 1 

E T l/ (X U:k )F 1 (X„ tk )J^(X„ jk t) = -p.v. / F 1 (x)J IM (xt)dx - -r lv [F 1 (z)J^(zt]} - 

k=i 1 Jo 1 

~T S {2-5 0vk )^ z= , k F 2 {z)J,{zt)^-±, (3.39) 
For example, it follows from here that for t < 1, Recr, Rez^ > — 1 

k=l v,k Z Jo Z Z L(<J + L) 



(for the value of integral see, e.g., [£(]]). For B = this result is given in 23]. In a similar 
manner taking fj, = v + m, 



.1 Ja(a\ z 2 + z?) 

v+m+l v V 1 



Fi(z) = 2 "+^^_v a>0 (3.41) 

(Z 2 +^) CT/2 

with Rei^ > and Rei^ + m > 0, from fl3.39f ) for a < 2 — t, Rea > Rev + m one finds 



^2T u (X^ k )X^ m+1 J u+m (X U:k t)—^- ' ^ a/2 - (3.42) 



( x 2 + z 2)-/2 a - (a 2 -i 2 ) (m+1)/2 



2 



00 „ +77l+1 . , Xja^g+g) f+i J m+1 (z iv ^^) 

* J »+rn{xt)—-^—— 7 ^dx = — (- Zl ) — , , 1W9 , a > t 



and the sum is zero when a < t. Here we have used the known value for Sonine integral |2Q| . 
If an addition to ( |3.37] ), fl3.38| ) the function F\ satisfies conditions 

Fi(a;e™ /2 ) = -e^Fiixe-™/ 2 ) (3.43) 

and 

\F 2 (z)\ < ei(x)e Clt \ y \ or \F 2 (z)\ < M\z\~ ai e tM , \z\ oo, (3.44) 



when the formula (5/7) (see below) with = may be applied to the integral on the right of 
(pi). This gives 0, |1 
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Corollary 1. Let F(z) be meromorphic function in the half-plane Rez > (except possibly at 
z = 0) with poles Zk, Rez^ > and ±iyk, Vk > \ u ,i)- If F(z) satisfy condition 

F(xe wi/2 ) = {-l) m+1 e uni F{xe- m/2 ) (3.45) 

with an integer m, and to one of inequalities 

\F(z)\ <si(x)e a M or \F(z)\ < M\z\~ ai e a °M, \z\ -» oo, (3.46) 

with a < min(t, 2 — t) = clq, £\{x) = o{x 1 / 2 ), x — > +oo, a.\ > 1/2, the following formula is valid 

oo 

^ T u {\ v ^k)F{\ v ^k)Jy+m{^v,kt) 
i; I 

= X ^ (2 " 5o ^ ReSz =^ { [ J v+™(zt)Y u (z) - Y u+m {zt)Uz)] £||} . (3.47) 



Recall that for the imaginary zeros X Ut k, m lhs of (3.47) the zeros with positive imaginary 
parts enter only. By using the formula ( 3.47 ) a number of Furier-Bessel and Dini series can be 
summarized (see, for instance, below). 



Remark. The formula (3.47) may be obtained also by considering the integral 



1 

7r Jc h 



H^Uzt)H^\z)-H^Uzt)H^{z) 



(i) 



F(z) 

J viz) 



dz, 



(3.48) 



where Ch is an rectangle with vertices ±ih, b ± ih, described in the positive sense (purely 
imaginary poles of F{z) / J v {z) and the origin are around by semicircles in the right half-plane 
with small radii). This integral is equal to the sum of residues over the poles within Ch (points 
Zk, ^u,k, (Rezk, ReA^fc > 0)). On the other hand it follows from ( p. 45 ) that integrals along the 
segments of the imaginary axes cancel each other. The sum of integrals along the conjugate 
semicircles give the sum of residues over purely imaginary poles in the upper half plane. The 
integrals along the remained three segments of Ch in accordance with ( |3.46 ) approach to zero 
in the limit b, h — ► oo. Equating these expressions for ( 3.4SQ one immediately obtains the result 

O- ■ 



From (3.47) for t = 1, F(z) = J u (zx), m = 1 one obtains |?], f20f 

OO 1/ 

E Tu(K,k)MK,kx)Ju+i(K,k) = \, < x < 1. (3.49) 

we obtain the Kneser- 



k=l 



By similar way choosing m = 0, F(z) = z,J u (zx)/ (z 2 — a 2 ) , B 



Sommerfeld expansion J20|: 

Jv{Xi/,kt)Ji/\^u,k x ) 



a 2 )J 2 +1 (K,k) 



tt J v (ax) 
4 J J a) 



[J u (at)Y v (a) - Y u (at)J u (a)] , < x < t < 1. (3.50) 



In ( 3.47 ) as a function F(z) one may choose, for example, the following functions 

-Mi/2 



■I )i 



(3.51) 
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for Ref < 22 -R- e W + n /2 + 2p + 3/2 — m — Sb ao , b < a , b = bf, 



1=1 



1=1 



Z' 



p-2n-\ 



A [i - Mkz)] 



1=1 



for Ref < 2n + 2p + 3/2 — m — 5b ao ; 
1 J! (:- ■ :f V ' 



( biy z 2 + z 2 ) , /x; > 0-half of an odd integer, 



(3.52) 



(3.53) 



for Ref < — /U/ + n/2 + 2p + 3/2 — m — <5fj ao ; 



i=i 



1 J\z^ [J fll+kl {aiz)Y n (biz) -Y^ l+kl (aiz)J n (biz)] , k t - integer, 



(3.54) 



i=i 



for Rev < n + 2p + 3/2 — m — /J |fc;| — <5g, ao , a = |a; — b;| < ap; 



i=i 



with p = v + m — 2p (p - integer), as well as any products between these functions and with 
Y\i (z 2 — of) Pl , provided the condition (|3.46j ) is satisfied. For example, the following formulae 
take place 



k=l J v+l{Jv,k) l=1 



W-2 



f_/'l_^TT_OL 



(r 



, o < t < 1, 



(3.55) 



71 



c = | a; — bi\ < t, at, b% > 0, Re/i/, > 0, Ref < n + 3/2 — 5, 
i=i 

Jp{3v,kt)Jv + l{\i V ,k) TT f1 r . _ A" +1 (1 - f 2 ") -A- 2 



ct ■ 



(3.56) 



X + J2 b l<t<l, A,6 Z >0; 
z=i 

( ji/,fc&) Jv+i (^3u,k) Jv (ju,kt) _ TrJ„ +1 (a\) JjApa) 

fc=l ° ) 3u,k J u+l\Ju,k) ^ I 

A + b < t < 1, A, b > 0, Re/i > -7/2 + <$ A+6)t , 



E 



[^(a) ^(at) - J^a^a*)] (3.57) 



where j v> k are zeros of J v (z). 

The examples of the series over zeros of Bessel functions we found in literature (see, e.g., 
[[7|, 20, 22]), when the corresponding sum was evaluated in finite terms, are particular cases 
of the formulae given in this section. 



4 Summation formulae over zeros of J u (z)Y v (\z) — Y v (z)J u (\z 

In this section we will consider the series over zeros of the function 

iAB , 



C (A, z) = J u (z)Y u (Xz) - Y u (z)J u (Xz), 



(4.1) 

where the bared quantities are defined as (|3.1|). Series of this type arise in calculations of the 
vacuum expectation values in confined regions with boundaries of spherical and cylindrical form. 
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To obtain a summation formula for these series let us substitute in (|2.11 



9{z) 



1 

27 



r( 2 )/ 



i^)J 



/(*) 



Hi X \z)Hl 2 \z) 



(4.2) 



where for definiteness we shall assume that A > 1. The sum and difference of these functions 
are 



s(z) - (-!)*/(*) 



fc = 1,2. 



(4.3) 



The conditions for GAPF written in terms of the function h(z) are as follows 

\h{z)\ < e 2 (x)e C2lyl or \h(z)\ < M\z\~ a2 e 2{x ~ lM , \z\ -> oo, z = x + iy (4.4) 



where c 2 < 2(A — 1), x 2<5so ^(^O - * for x — > +oo, a 2 > 25bo- Let 7^ be zeros for the 
function C^ B {\, z) in the right half-plane. In this section we will assume values of v, A, B for 
which all these zeros are real and simple, and the function H„(z) {&v [z)~) has no zeros in 
the right half of the upper (lower) half-plane. As we will see later these conditions are satisfied 
in physical problems considered below. For real v, A, B the zeros 7^ are simple. To see this 
note that the function J u (tz)Y u (z) — Y u (tz)J u (z) is cylinder function with respect to t. Using 
the standard result for indefinite integrals containing any two cylinder functions (see p0| , |2l[| ) it 
can be seen that 

' X t [J u (tz)Y u (z) - Y v {tz)J u {z)} 2 dt 



1 

where we have introduced the notation 



ir 2 zT^ B (X,z) 



T u AB (X,z) 



On the other hand 



fJM. 

\j*(Az) 



A 2 + B 2 {X 2 z 2 -v 2 ) 



A 2 



B 2 {z 2 -v 2 ) 



9 ^AB 



C™{X,z) 



7u,k- 



dz~ u v "'~ / vrT^(A,z) J v {z) ' 

Combining the last two results we deduce that for real u, A, B the derivative (|4 
and hence the zeros z = 7^ k are simple. By using this it can be seen that 



Res 



7T 



2 =7, 



(4.5) 

(4.6) 

(4.7) 
is nonzero 

(4.8) 



Hence if the function h(z) is analytic in the half-plane Rez > a > except at the poles z^ 
lv,i) and satisfy to the one of two conditions (4.4), the following formula takes place 



\ £ T AB (X, lv>k )h( lv , k ) + r 2u [h(z)] - p.v. / 
2 k=n Ja 



h{x)dx 



a+ioo 



Rj l \Xz) h(z) 
Hi l \z) C? B (X„ 



-dz 



H {2) {Xz) 



J 2 (x)+Y 2 (kx) 

Hz) 



Si 2 \z) C**{\z) 



dz. 



(4.9) 



Here we assumed that the integral on the left exists, 7^-1 < a < 7^ n , 7^ m < b < ^y u ,m+i, 
a < Rezk < b, Rezk < Rezk+i, and the following notation is introduced 



r 2v [h(z)] 



7r^Resi m2fc= o 
k 

+7r ^ Res^ 

fc,/=l,2 



J u (z)J u (Xz) + Y u (z)Y u (Xz) h(z) 

P u {z)+Y 2 {z) C^ B (X,z) 

Sl l) {Xz) h(z) 



+ 



) ; lmz fc <0 



H, 



(0, 



C^ B (X,z) 



(4.10) 
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The general formula (4.9) is a direct consequence of GAPF and will be as starting point for the 
further applications in this section. In the limit a — ► one has [11, 12]: 

Corollary 2. Let h(z) be analytic function for Rez > except the poles Jui), Rez^ > 

(with possible branch point z = 0). If it satisfies one of two conditions (U-4\) an d 



and the integral 



h(ze wi ) = -h{z)+o{z- 1 ), 
h(x)dx 

""J a J 2 (x)+Y 2 ( X ) 



0, 



p.v. 



(4.11) 
(4.12) 



exists, then 



hm \ Y, h (7u,k)T* B (\, lu , k )+r 3u [h(z)]-p.v. / 



h{x)dx 



7T 



Res 



2 = 



K v {\x) 



J 2 (x) + Y 2 (x) $ 



dx 



h(z)Si 1 \\z) 

C^ B (X,z)Hl 1 \z)\ 4 ^o K v (x) K v (x)I v {\x) - K v (\x%(x) 



(4.13) 



In the following we shall use this formula to derive the regularized vacuum energy momentum- 
tensor for the region between two spherical and cylindrical surfaces. Note that ( |4.13| ) may be 
generalized for the functions h(z) with purely imaginary poles ±iyk, Uk > satisfying condition 



h(ze T ' 



-h(z) + o ((z =f 



Wk) 



±iyk- 



(4.14) 



The corresponding formula is obtained from ( 4.13Q by adding residue terms for z = iyk in the 
form of ( [4,16| ) (see below) and taking the principal value of the integral on the right. The 
arguments here are similar to those for Remark after Theorem 3. 
By the way similar to ( |3.18| ) one has another result [11, |l2^ : 

Corollary 3. Let h{z) be meromorphic function in the half-plane Rez > (with exception the 
possible branch point z = 0) with poles z^, ±iyk (j^ jvi), Rezk,yk > 0. If this function satisfy 
condition 

h{xe™l 2 ) = -h{xe~ m/2 ) 



(4.15) 



and the integral ( \4-lg ) exists then 



f 7r 2 m rb 

hm \—J2 h (l »,k)Tf B (A, 1Ujk )+r 2u [h(z)]-p.v. \ 
~+oo { 2 ^ Jo 



h(x)dx 



J 2 (x)+Y 2 (x) 



Vk=o,iyk 



°0r, k 



Res 



z=Vk 



H$\\z) 



h(z) 



SP(z) 



C^(A, 



(4.16) 



where in the Ihs j u ^ m < b < ^u^+i- 



Let us consider a special applications of the formula ( |4.16| ) for A = 1, B = 0. The general- 
izations of these results for general A, B under the conditions given above are straightforward. 

Theorem 4. Let the function F(z) be meromorphic in the right half-plane Rez > (with the 
possible exception z = 0) with poles z k , ±iyk (£ lut), yk, Rez^ > 0. // it satisfy condition 



F(xe ni / 2 ) = (-l) m+1 F(xe 



m+l 



-iri/2\ 



(4.17) 
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with an integer m, and to the one of two inequalities 

\F(z)\ < e(x)e aiM or \F{z)\ < M\z\~ a e a ^ y \ \z\^oo, (4.18) 
with a\ < min(2A — a — 1, cr — 1) = cr > 0, e(x) — > /or x — ► +00, q > 1, then 

1 ro x \p y„(Az)J„ +m (gz) - J^Az^+^gz) 



Proof. As a function h{z) in ( 4.16| ) let us choose 



= F(z) [J v {z)Y v+m {az) - Y u {z)J v+m {az)) , (4.20) 
which in virtue of ( 4.18 ) satisfy condition (|4.4j ). The condition ( ^4.15|) is satisfied as well. Hence 



h(z) satisfy conditions for Corollary 3. The corresponding integral in (4.16) with h(z) from (4.20) 
can be calculated by using the formula ( [7.7| ) (see below). Putting the value of this integral into 
( |4.1GD after some manipulations we receive to fl4.19|) , ■ 



Remark. The formula ( 4.19| ) may be derived also by applying to the contour integral 



f Y v (Xz)J v+m ((Tz) - J u {Xz)Y u+m {crz) 

J Ch Ju(z)Y u (Xz) - J v {Xz)Y v {z) * [Z)dZ l4 - 2ij 

the residue theorem, where Ch is a rectangle with vertices ±ih, b ± ih. Here the proof is similar 
to that for Remark to the Corollary 1. ■ 



Formula similar to ( 4.1 9| ) can be obtained also for the series of type J2k=iG(lv,k)Jfj,(7u,kt) 



by using (|4.16[) 



As a function F{z) in (4.19) one can choose, for example, 



• function fl3.51| ) for p = m — 2p, J2i h < a 2, m < 2p + J2i Re//; + n/2 + 1, p - integer; 

• function ( |3.52| ) for p = m — 2p, bi < 02, m < 2p + 2n + 1; 

• function ( |3.54| ) for p = m — 2p, a/ > 0, Re/i; > (for Re/i/ < 0, ki > |Re/^|), 
YZ=i \ai -b t \ < a 2 , m <2p + n- J2i \k\ + 1. 

For F{z) = 1/z, m = one obtains 

h — — = ~ a--i ' (422) 

By similar way it can be seen that 

ll,k [Jy{lu,k)Yu{o-lv,k) - Y u (^ k )J u (a^ k )] = 1 Y v {\c)J v (ac) - J v (Xc)Y v (crc) 
^ (ll, k -c 2 )[JK7u,k)/JK^k)-l] _7r Jv{c)Y v {Xc) - J v {Xc)Y v {c) ' ' 



k=i 

00 



E -Jy \ lu,kjiy\v lv,k) ~ *-v\ lu,k)-Jy\v lu,k) rr -ft r /. \ _ 
J 2 Aly,k)/J^,k) ~ 1 l\ % > k J " l[bnU ' k) ~ 
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where Rec > 0, b x > 0, A > a > 1, ft ^ — 1, -2, . . .. 

So far in this section we have considered series over zeros of the function C^ B (X, z). The 
similar results can be obtained also for the series containing zeros of the function 

C lv {\, z) = j' v {z)Y v {Xz) - Yl(z)J v (\z), (4.25) 

(on properties of zeroes of these function see [p0| , [j], ^]). The corresponding furmulae for the 
zeros jijy^k of this function can be obtained from those for C^ B (X, z) by replacements 

T* B (A, Z) - j2 (z)/J 2 (Az) _ 1 + z ,2 /z 2 ( 4 - 26 ) 



/(*)-/'(*), f(Xz)^f(Xz), f = J u ,Y u ,Hj ) 1 ' 2 \l u ,K u , Cf^C 



The physiacl apllications of the formulae derived in this section will be considered below in 
Sections 9 and 12. 



5 Applications to integrals involving Bessel functions 

The applications of GAPF to infinite integrals involving some combinations of Bessel functions 
lead to the interesting results 11, 12]. First of all one can express integrals over Bessel functions 
through the integrals involving modified functions. Let us substitute in the formula ( 3.1 8| ) 



f(z) = F{z)J v {z). 



(5.1) 



For the function F(z) having no poles at z = X V) t the sum over zeros of J u (z) is zero. The 
conditions (|J) and (|t|) may be written in terms of F{z) as 



\F{z)\ < ei(x)e Cl|s/| or \F(z)\ < M\z\~ ai e M , \z\ -»• oo, 
with ci < 1, x 1 / 2 ~ Sbo Ei(x) — > for x — > oo, a\ > ao = 3/2 — 5bo, an d 

F{ze Ti ) = -e™F{z) + o ( z l R H-i x 



(5.2) 



(5.3) 



Hence for the function F(z) satisfying conditions (|5.2| ) and ( |5,3| ) it follows from ( 3.1S| ) that 



7T 



p.v. / F{x)J v {x)dx = r lv [F{z)J v {z)] + -Res z=0 F(z)Y u (z) + 



+ - 



vr Jo 



K u (x) \e- uwi/2 F(xe ni/2 ) + e U7ri/2 F(xe- wi/2 ) 



dx. 



(5.4) 



In expression ( |3.S| ) for r\ v the points Zk are poles of the meromorpic function F(z), Rezk > 0. 
On the base of Remark after Theorem 3 the formula ( |5.4[ ) may be generalized for the functions 
F(z) with purely imaginary poles Hyk, Vk > and satisfying condition 



F(ze^) = -e» m F(z) +o((z^ iy k )- 1 



(5.5) 



The corresponding formula is obtained from ( |5.4| ) by adding residue terms for z = iyk in the 
form of ( |5.7D (see below) and taking the principal value of the integral on the right. 



The same substitution (5.1) with the function F{z) satisfying the conditions (|5.2j) and 



F(xe m/2 ) 



^(xe 



-7ri/2\ 



(5.6) 



for real x, into the formula (3.33) yields the following result 
p.v. 



F{x)J v {x)dx = r lv [F{z)J v {z)]+- £ (2 - 5 0r]k ) Res z=Vk F(z)H^ (z). (5.7) 

Vk=0,iy k 
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In ( p.9| ) now summation is over the poles Zk, Rez/t > of the meromorphic function F(z), and 
=t*2/fcj Vk > are purely imaginary poles of this function. Recall that the possible real poles of 
F(z) are such, that integral on the left of fl5.7| ) exists. 

For the functions F(z) = z u+1 F(z), with F(z) being analytic in the right half-plane and 
even along the imaginary axis, F(ix) = F(—ix), one obtains 



x u+l F(x) J v {x)dx = 0. 



(5.8) 



This result for B = (see ( |3.1|) ) have been given previously in [24]. The another result of [24] 
is obtained from (5.7) choosing F{z) = z u+l F(z)/(z 2 — a 2 ). 



Formulae similar to ( |5.4| ) and fl5,7| ) can be derived for Neumann function Y u (z). Let for the 
function F(z) the integtral p. v. J °° F(x)Y u (x)dx exists. Let us substitute in the formula ( |2.11| ) 



f(z) = Y u (z)F(z), g{z) = -iJ u (z)F(z) 



(5.9) 



and consider the limit a — > +0. The summands containing residues may be presented in the 
form 



R[f(z),g(z)} = 7T £ Res lmZk>Q H^ (z)F(z) + n £ Resi mZfc<0 ^ 2) (z)F(z) + 

k k 

+ 7r^Resi mZfc=0 J ! ,(2;)i ? (2;) = r 3u [F(z)\, 

k 



(5.10) 



where Zk (Rez^ > 0) are the poles of F{z) in the right half-plane. Now the following results can 
be prooved by using ( 2.11 ): 



condition ( 5A ) then 



1) If the meromorphic function F(z) has no poles on the imaginary axis and satisfy the 

K v {x) e- vm/2 F{xe ni/2 ) - e um/2 F(xe~ ni/2 ) 



p.v. / F(x)Y u (x)dx = r 3iy [F(z)} 



vr Jo 



d,x 
(5.11) 



and 

2) If the meromorphic function F(z) satisfy the conditions 

F(xe m l 2 ) = e um F{xe- m / 2 ) 

and ( \5.Q ) then one has 

/•OO 

p.v. / F{x)Y u {x)dx = r 3u [F(z)]+Try2Re Sz=iyk Hi 1 \z)F(z 
Jo k 

where Hyk, Vk > are purely imaginary poles of F(z). 



(5.12) 



(5.13) 



From ( |5.13D it directly follows that for F(z) = z u F(z), with F(z) being even along the 
imaginary axis, F{ix) = F(—ix), and analytic in the right half-plane 



p.v. / x"F(x)Y v (x)dx = 0, 



(5.14) 



if the condition (|5.2|) takes place. 
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Let us consider more general case. Let the function F(z) satisfy the condition 

F(ze- ni ) = -e- Xni F(z) (5.15) 
for argz = n/2. In GAPF as functions f{z) and g{z) we choose 
f(z) = F(z) [J v {z) cos 5 + Y v (z)sm 5] 

g(z) = -iF{z) [J u {z) sin 5 - Y u (z) cos 5, 6 = (A - u)ir/2] , (5.16) 

with g(z) — (—l) k f(z) = H^\z)F(z) exp[(— l) k i6], k = 1, 2. It can be seen that for such a choice 
the integral on rhs of (|2.11| ) for a — ► is equal to 

ni J2 Res z=Vk Hi 1 \z)F(z)e- iS , (5.17) 



where Hyk, Uk > 0, as above, are purely imaginary poles of F(z). Substituting ( 5.16j ) into ( 2.11 ) 
and using ( ^Oj ) we obtain [11, 12 1 

Corollary 4. Let F{z) be meromorphic function for Hez > (except possibly at z = 0) with 
poles Zf., ±iyk! yfcjRezjt > 0. If this function satisfies conditions \5.t\ ) (for B = 0) and f |5. 
then 

poo I 

p.v. / F(x)[J u (x)cos5 + Y u (x)sinS]dx = Tri<y^Resi mZk>0 Hl 1 \z)F(z)e' i5 - 
Jo \X 

- Resi mZk<0 H^ (z)F(z)e iS - i ^ Resi m2fc=0 [ J„(z) sin 5 - Y u (z) cos 5] F(z) + 



+ £ Re Sz=Vk H^(z)F(z)e- i5 \, (5.18) 

Vk=iVk ) 

where it is assumed that integral on the left exists. 

In particular for S = irn, n = 0, 1, 2... the formula ( |5.7D follows from here in the case B = 0. 
One will find a great many particular cases of the formulae (|5.7|) and (|5.18j) looking at the 



standard books and tables of known integrals with Bessel functions (see, e.g., [0, £C], 22, 25 



27, Eq, E9|). Some special examples are given in the next section. 



6 Integrals involving Bessel functions: Illustrations of general 
formulae 

To illustrate the applications of the general formulae from previous section first of all consider 
integrals involving the function J u (z). Let us introduce the functional 

fOO 



poo 

Ap, m [G(z)] = p.v. / z v - 2m - l G{z)J u {z)dz (6.1) 
Jo 



with m being an integer. Let F\{z) be an analytic function in the right half-plane satisfying 
condition 

F^xe™' 2 ) = F^xe-^ 2 ), Fi(0) ^ (6.2) 
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(the case when F\(z) ~ z g , z —* with an integer q can be reduced to this one by redefinitions 



of F\(z) and m). From ( |5.7| ) the following results can be obtained [11, 12 1 

2m 

.dz J 



A, 



A um [Fi(z)\ 
Fi{z) 



2 = 



z z — a z 



7T 



v-2m-2 



F v (a)Fi(a) + AW 



z 2 — a 2 



,j/-2m-4 r 



r y,(a)Fi(a) 



-l) m ^(a)Fi(ia) 



+ 



+4(0) 



Fi{z) 



A, 
A„ 



[(z 2 -c 2 ) p+1 
Fi(z) 
z 2 + a 2 ) p+1 



(-l) m +P +1 / d y 
2P-p\ Voda/ 



F 1 (c) J F/( 1 )(c)l+4 ,) 



Fi(z) 



f z 2 _ C 2)P+1 



(6.3) 
(6.4) 

(6.5) 
(6.6) 



v-2m-2 



^(a)Fi(ae 



+ 



+ 4(0) 



Fi(z) 



(z 2 +a 2 ) 



2\P+ 1 



(6.7) 



and etc. (note that A U m = for m < 0). Here a > 0, < argc < 7r/2, and we have assumed that 
Ref > 0. To secure convergence at the origin the condition Hev > m is necessary. In the last 
two formulae we have used the identity 



d \ p 
~dz) 



zF(z) 

(z + b)P+ 1 



1 



2PH 



d V 



(6i 



Note that ( |6.7| ) can be obtained also from ( |6 .6|) in the limit Rec — ► 0. For the case F% = 1, m 

[]) as a function -Fi(z) we can choose: 



1 of (|6.7D see, for example, [ 20 1 . In (|67 
• function (|3.51| ) for p = 1, Rez/ < £ Rep^ + 2m + (n + 1)/2 - <5&i + Sbo, b = £6; < 1, k > 0; 



• function ( gjjg) with p = 1, Rez; < 2(m + n) + 1/2 - <5 61 + <5 B o, & = E k < 1; 

• function fl3.53|) for p = 1, Rev < 2m — J2 Rep; + (n + l)/2 — 5bi + #bo 5 pz > is half of an 
odd integer, b = < 1; 



• function ( gjg ) for Rezv < 2m + n - £ + 1/2 + 
ki - integer. 



hi, a = J2\a>i 



<l,a t > 0, 



Here we have written the conditions for (|Q|), The corresponding ones for ( |6.4| ), fl6.5| ), ( |6.6| ) , (p/7]) 
are obtained from these by adding on the rhs of inequalities for Rez/, respectively 2, 4, 2(p + 1), 
2(p + 1). In ( |6.3| )-( |6?7| ) we can choose also any combinations of the functions (3.51 )-( |3.54| ) with 
appropriate conditions. 

For concrete evaluations of A^m in special cases it is useful the following formula 



d 



2/d 



lim 

z^o \dz J 



h{z) = (2m-l)!!lim ( -- 

z^o \zdz 



fi(z), 



valid for the function f±(z) satisfying condition fi(—z) = f\(z) + o(z 



2m \ 



(6.9) 

0. From here, for 



instance, it follows that for z — ► 

d x 2m 



-J [z^Y u (bz)F 1 (z)} 



P— i)"a^£( t ja t sb H is) (6 ' 10) 
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where we have used the standard formula for the derivative {d/zdz) n of cylinder functions (see 
Q). From (O) one obtains (B = 0) 



v v-2m-l 



7T 



1=1 

d 



2 , 2 

Z +Z l 



-Mi/2 



2 m + 1 m! \zdz 



J m( b iyz 2 + zf)dz = 

^)n(z a +^)" w/2 j tt (6 I v^+?) 



(6.11) 



2=0 



for to > and the integral is zero for to < 0. Here Rev > 0, b = J2i^i — 1> h > 0, 
to < Rev < Yl?=i R- e ^ + 2to + (n + 3)/2 — <5&i. In particular case to = the Gegenbauer 
integral follows from here [7J. In the limit Z[ — > from ( 6.1 1| ) the value of integral 
j°c z u-2m-ij^ g-mj^dz can be obtained. 

By using ( |5. 18 ) the formulae similar to fl6.3| )-(|67^) may be derived for the integrals of type 

B u [G(z)]=p.v. G(z) [J u (z) cos 5 + Y v (z) sin 5\dz, S = (A - i/)tt/2. (6.12) 
Jo 

It directly follows from Corollary 4 that for function F(z) analytic for Rez > and satisfying 
conditions (5.2) and ( 5.15| ) the following formulae take place 

(6.13) 
vrF(a) [J u {a) sin 5 - Y v {a) cos S] /2, (6.14) 

^F(a) [J u (a)B3n6-Y v (a) cos 5] + ^K u (a)F(ia)e~ iXn/2 ,(6.15) 





B,[F(z)] = 


B„ 


" F(z) - 




z 2 — a 2 




B u 


f : 




U 4 -a 4 _ 




Fi{z) 1 




> 2 - 


-c 2 ) p+1 _ 




Fi(z) 




_( z 2 +a 2)P+l_ 





TTl 



P r 



2P+1 . p \ \cdcj 

(-i)p +1 f d y 

2P -p\ \adaj 



'c- x F{c)H(% 



ar l F{ae™l 2 )K v {a) 



(6.16) 
(6.17) 



where a > 0, < argc < tt/2. To obtain the last two formulae we have used the identity 
fl6.8D - The formula ( |6.13| ) generalizes the result of |£4j (the cases A = v and A = v + 1). From 
the last formula taking F{z) = we obtain result given in |£Q]. In ( |6.13| ) - (|6.17| ) as a 

function F(z) one can choose (the constraints on parameters are written for the formula J6.13 ); 
the corresponding constraints for (|6.14| ), ( |6.15| ), ( |6.16| ), ( |6.17 ) are obtained from given ones by 
adding the summands 2, 4, 2(p + 1), 2{p + 1) to the rhs of inequalities, correspondingly): 

• function ( |3.5lD for p = A, \Rev\ < Re/? < Repi + (re + 3)/2 — S^i, b = J2i h < 1; 

• function ( |3.52 ) for p = A, |Rez/| < Rep < 3/2 — Sbi, b = J2i bi < 1; 

• function 

n n 

zP 1 II [ J ^i+kii.aiz)Y N (biz) - Y )ll+kl (aiz)J N (biz)} , A = p + ^^, a t > 0,(6.18) 
1=1 i=i 

\Rev\ + ^2\h\ < Rep < n + 3/2 - 5 cl , c = ^ \ai - b t \ < 1, Re/i; > 

(for Repi < one has k\ > \Repi\). 

Any combination of these functions with appropriate conditions on parameters can be choosed 
as well. 
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Now consider integrals which can be expressed via series by using 

Q and fU8P . In Q 

let us choose the function 

F * = , u ; , 6.19 
sinn irz 

where F\{z) is the same as in the formulae (|6.3| ) - (|6.6|), As the points ±i, ±2i, . . . are simple 
poles for F(.z) from (5.7) one obtains 



co z v—2m 

o sinh (7rz) 



F 1 {z)J v {z)dz = A^l 



zF x {z) 



sinh (Vz) 



+ - Yi-ir+^v-^K^F^ik), (6.20) 
, — i 



k=l 



where Abm[f{z)] is defined by fl6.3|) and Ref > m. The corresponding costraints on F\{z) follow 
directly from ( |5.2| ). The particular case of this formula when F\(z) = smh(az) / z and m = — 1 is 
given in [2C[ ]. As a function -Fi(z) here one can choose any of functions ( p.51| )-( |3J>4] ) with p = 1 
and a, J2i bi < 1. From ( 6.20| ) it follows that 



oo 2m 



sinh (7rz 



z=i 



+ 



+ _^(_l)m+fe iCi/(A , ) jjAj-wj^^), 6, > 0, tt - £&, > 0, Rei/ > m. (6.21) 



fc=i 



In similar way from ( 5.1q ) it can be derived the following formula 



zF(z) 
sinh (irz) 



2i 



JJ z) cos d + YJz) sin S\dz = —e- iXn/2 Y{-l) k kK u {k)F(ik). (6.22) 

7T 



k=\ 



Constraints on the function F(z) immediately follow from Corollary 4. Instead of this function 
we can choose the functions (jjTftiD , ( gjp , (jjf ). 

As it have been mentioned above adding residue terms -KiKes z= i yk F{z)H^' \z) to the rhs 
of (|5.4| ) this formula may be generalized for the functions having purely imaginary poles Hyk, 
Uk > 0, provided the condition (|5.5|) is satisfied. As an application let us choose 



F(z) 



z v F x (z) 

2irz/b _ i ' 



F 1 (-z) = F 1 (z), b>0 



(6.23) 



with an analytic function F\(z). The function ( |6.23 ) satisfy condition ( |5.5| ) and have poles ±ikb, 
k = 0, 1, 2 . . .. The additional constraint directly follows from (|5.2| ). Then one obtains 



o e 



V J ( \ q OO 

^-£^L Fl {x)dx = -Y' {bkY K v {bk)F^bk) 



k=0 



1 



x u K u (x)F 1 (ix)dx, (6.24) 



where the prime indicates that the m = term is to be halved. For the particular case F\(z) = 1, 
using the relation 



£ '{bkYK v {bk) = ^-<FY(y + 1/2) £ ' 



-i -I/-1/2 



+ 1 



(6.25) 



and the known value for the integral on the right, we immediately obtain the result given in 



(6.26) 



The relation ( |6.25| ) can be proved by using the formulae 
2T(i/ + l/2) f°° cos ztdt 



K v {z) 



sftz v Jo (t 2 + l)H-i/2 ' 



J2 e ikz = 2n 2 <5(*-27rn), 



fc=— oo 



(for the integral representation of Macdonald's function see [20]), S(z) is the Dirac delta function. 
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7 Formulae for integrals involving J v (z)YJ\z) — Y v (z)JJXz) 

In this section we shall consider the applications of GAPF to the integrals involving the function 
J u (z)Y fl (Xz) — Y v (z) J ]i(Xz) . In the formula ( p. 11 ) we substitute 



1=1 



H®(z) 



2i 



(7.1) 



For definiteness we consider the case A > 1 (for A < 1 the expression for g{z) have to be 
choosen with opposite sign). The conditions (|2j]) and fl2~lc| ) are satisfied if the function F(z) 
is constrained by the one of the following two inequalities 



\F(z)\ < e(x)e c ^, c < A - 1, e(x) -> 0, x 



-oo 



(7.2) 



or 



\F(z)\ < M\z\- a e {x ~ 1)lyl , a > 1, \z\ -> oo, z = x + iy. (7.3) 
Then from ( |2.1l| ) it follows that for the function F(z) meromorphic in Rez > a > one has 

J u (x)Y^(Xx) - J^{Xx)Y v {x) 



p.v. 
1 



Jd{x)+Y^ 



X 



+ 



2i 



F{z) H ^} XZ) dz 



H<P(z) 



-F(x)dx = r llu/ [F(z)] + 



H^{z) 



(7.4) 



where we have introduced the notation 

2 



ri^[F(z)] = -J2^l 



Im2:j.=0 



WE 



HS\\z) 



1=1 Hl l) (z) 



) ; Im2; fe <0 



k 1=1 



F(z 



H$\\z) 
H«\z) 



(7-5) 

The most important case for the applications is the limit a — > 0. The following statements take 
place UJ, |TJ: 

Theorem 5. Let the function F{z) be meromorphic for Rez > (except the possible branch 
point z = 0) with poles z^, ±iyk (yki^Zk > 0). If this function satisfy conditions ( \7.3( ) or ( \7.3{ ) 
and 



Fixe 71 ' 1 / 2 } 



F(xe 



-m/2\ 



(7.6) 



then for values of v for which the function H^\z) (H^\z) ) have no zeros for < argz < tt/2 
(—it/2 < argz < 0) the following formula is valid 



p.v 



J^Y^Xx) -Y^J^Xx) 



JZ(x)+Y?(x) 
+1 E (2-5o Vk )Re Sz=Vk F(z)^^, 



Vk=0,iVk 

where it is assumed that the integral on the left exists. 

Proof. From the condition ( [7.6|) it follows that for argz = ir/2 

F(zt), zi = 



F(x)dx = r lflu [F(z)} + 



(7.7) 



H^(Xz) H^{X Zl ) 



H£\z) 



-F z) 



h1 2 \ Zi ) 



(7.8) 
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and that the possible purely imaginary poles of F{z) are conjugate: ±iyk, Dk > 0. Hence in rhs 
of ( |7.4p in the limit a — > the term in the square brackets may be presented in the form (it can 
be seen similarly to ( 3,32| )) 



F(z)dz (7.9) 



with the same notations as in (|3,32j ). By using (|7.8| ) and the condition that the integral converges 
at the origin we obtain 

I Ht^M F{z)dz+ i E¥^M F{z)dz = {2 _ s 071k)mRes ^££) F(z)i 

(7.10) 

where fip (0) = xj, Vt±{iy^) = C p (±iyk). By using this relation from ( |7.4[ ) we receive the 
formula (^). ■ 

Note that one can write the residue at z = in the form 

Hp\\z) J u (z)J,(\z)+Y v (z)Y,(\z) 
Z= °ltf\z) {) = 2=0 ■£(*)+*?(*) F{Z) (7 - 11} 



as well. Integrals of type fl7.7p we have been able to find in literature (see, e.g., |7j |23|, |25|, |27| ) 
are special cases of this formula. For example, taking F(z) = J u (z)Y u+ \(\' z) — Y v {z) J u+ i(X'z) 
for the integral on the left in Q we obtain -X^X'-"- 1 for A' < A and A^A'^" 1 - X^X'-"' 1 
for A' > A pjf|. By taking z 2m+1 /(z 2 + a 2 ), z 2m+1 /(z 2 - (?) as F\[z) for jj, = v and integer 
m > one receive 

J v {x)Y v {Xx) - Y u {x)J v {\x) x 2m+1 _ 2m TT K u (Xa) 

J 2 {x)+Y 2 (x) x 2 + a 2dX ~ [ L) a 2 K>) 1 ] 

[°° J u {x)Y v {Xx)-Y u {x)J v {Xx) x 2m+1 7T 2m J,(c)J,(A C )+y,(c)n(A C ) 

p ' v 'io J2(x) + y, 2 (x) * 2 -c 2Cte 2 C J 2 ( c ) + y 2 ( c ) u ' ^ 

where c > 0, A > 1. The particular cases of this formula for v = m = are given in J25jj . In 
( |7 . 1 2| ) taking the limit a — > and choosing ra = 0we obtain the integral of this type given in [23] . 
In ( |7.7| ) as a function F(z) we can choose ( |3.51| ), ( p. 52 ), ( 3.54| ) (the corresponding conditions 
for paremeters directly follow from (7J2) or ( |7.3| )) with p = \i — v — 2m (m - integer), as well as 
any products between them and with Il^=i( z2 — c f)~ kl - For instance, 



r Mx)Y u (Xx) - Y u (x)J u (Xx) " J* ibi^x 2 + zf) dx n « 

Jo idw^m l\ (x 2 +z?r/ 2 x~^L[ J » AbiziU1AA) 

n n 

b h Rev > 0, Rezi > 0, A > 1, ^ Re^ + n/2 + 1 > 5 biX -i, b = ^h < A - 1 

l=i i=i 

As another consequence of ( |7.4[ ) one has: 

Theorem 6. Let F{z) be meromorphic in the right half-plane (with possible exception z = 0) 
with poles Zk, Rez^ > 0, and satisfy conditions n?J\) or (\7.Q) and 



F{ze™) = -e^- u)7Ti F(z) + o f z M-\^\~A , z _» o, (7.15) 
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then for values of v for which the function H^\z) (H^\z) ) have no zeros for < argz < tt/2 
(—it/2 < argz < 0) the following formula takes place 

f°° JJx)YJXx)-YJx)JJXx) „, , , r „, tt HjpfXz)^,. 

P.v. I {) % { J )+Y ^ A > F{x)*c = V[FW] + -Jto^-J^-lFiz) + 



1 f°°K„(\x) 



2 Jo K v {x) 
provided the integral on the left exists. 



(^W2 F ( xe W2) + e (n-v)m/2 F ( xe -vi/2j\ ^ A > 1, (7.16) 



Proof. This result immediately follows from ( |7,4|) in the limit a — > and from ( 7.10 ) with 

m = o. ■ 

For example, by using ( |7,16| ) one obtains 

n n n 

^Refii + \Reu\ > |Re//| - 1, 6 = ^ < A - 1, > 0, n > 4,A-i, Ms = f - A 4 + 

Such relations are convenient in numerical calculations of integrals on the left as the subintegrand 
on the right at infinity goes to zero exponentially fast. 

We have considered the formulas containing J u (z)Y ll (Xz) — Y u {z)J^{Xz). The similar re- 
sults can be obtained for integrals containing the functions J' u (z)Y fl (\z) — Yl(z)J^(Xz) and 
J' u (z)Y^Xz)-Yl(z)J^Xz). 

8 Applications to the Casimir effect. Vacuum energy density 
and stress inside a perfectly conducting spherical shell 

In this and next sections we shall consider applications of the summations formulae obtained 
in previous sections to the physical problem, namely the Casimir effect. In what follows on the 
example of spherical and cylindrical geometries we will show that the using of GAPF allows 
to obtain the regularized values of physical quantities in cases then the explicit dependence of 
eigenmodes on quantum numbers is complicated and irregular. 

Historically the investigation of the Casimir effect for a perfectly conducting spherical shell 
was motivated by Casimir semiclassical model of an electron. In this model Casimir suggested 
that Poincare stress to stabilize the charged particle could arise from vacuum quantum fluctua- 
tions and the fine structure constant can be determined by a balance between the Casimir force 
(assumed attractive) and the Coulomb repulsion. However, as it have been shown by Boyer 
p0|] , the Casimir energy for the sphere is positive, implying a repulsive force. This result has 
later been reconsidered by a number of authors |3l], |32|, |33|]. More recently new methods have 



been developed for this problem including a direct mode summation [38, [39[ and zeta function 
[ ^0| , 41, 42 1 approaches. However the main part of studies have focused on global quantities 
such as total energy. The investigation of the energy distribution inside a perfectly reflecting 



spherical shell was made in [34| in the case of QED and in for QCD. The distribution of 
the other components for the electromagnetic EMT inside as well as outside the shell can be 
obtained from the results of p6|, |37[] . In these papers the consideration was carried out in terms 



of Schwinger's source theory. In |l3|, 14, |l^] the calculations of the regularized values for vev 



of the EMT components inside and outside the perfectly conducting spherical shell are based 
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on the generalized Abel-Plana summation formula. Our consideration below is based on this 
approach. 

The main quantities we will consider here are vacuum expectation values (vev) of the energy- 
momentum tensor (EMT) for the electromagnetic field inside a perfectly conducting spherical 
shell of radius a. It may be obtained by using the standard formula of mode summation 0. 53) 



(o|r«|o) = {«£->(*), *«(*)} 



where bilinear form g} for a field $ is given by the classical EMT. Here |0) is the amplitude 

of the vacuum state, j\Pa (^)| is a complete set of the positive and negative frequency solutions 
to the field equations, satisfying the boundary conditions, and subscript a may contain discrete 
and continous components. 

In the case of the electromagnetic field inside the perfectly conducting sphere, by using 
Coulomb gauge for vector potential A, the corresponding system of solutions, regular at r = 0, 
can be presented in the form 

a -i R ( ^ / ji(ur)X* m e-*" if A = 

Aa = lV ^'^(^VxfeHXtele-** if A = l ' « = M™ A )> ( 8 - 2 ) 

where A = and 1 correspond to the spherical waves of electric and magnetic type (TM and 
TE - modes). They describe photon with definite values of total momentum I, its projection m, 
energy u> and parity (— l)'+ A + 1 (units h = c = 1 are used). Here the vector spherical harmonics 
have the form 

X, m (0, <p) = -i-^=^=Y lm (9, <p), 1^0, (8.3) 



with Y[ m being spherical function, and ji(x) = v /7r/2xJ; +1/ / 2 (x) is spherical Bessel function. The 
coefficients 0i(a,u) are determined by the normalization condition 

2ir 

dVA a -A* a , = —5 Qa ,, (8.4) 

where the integration goes over the region inside the sphere. Using the standard formulae for 
vector spherical harmonics and spherical Bessel functions (see, for example, p3) one finds 



/3f(a,uj) = 8uj 3 T u (ua)/a, v = 1 + 1/2, (8.5) 



where T u (z) is defined in Q3.8[ ). 

Inside the perfectly conducting sphere the photon energy levels are quantized by standard 
boundary conditions: 

rxE = 0, r-B = 0, r = a, (8.6) 

where E and B are the electric and magnetic fields. They lead to the following eigenvalue 
equations with respect to u 

3l{ur)\ r=a = if A = (8.7) 

JrfalMWa = if A = l (8.8) 

It is well known that these equations have infinite number of real simple roots |2l|] . 

By substituting the eigenfunctions into (8.1) with the standard expression of the electromag- 
netic EMT and after the summation over m by using the standard formulae for vector spherical 
harmonics (see, for example, |j44}] ) one obtains 

(0|T ifc |0) = diag (e, -p, -p ± , - P± ) (8.9) 
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(here index values 1,2,3 correspond to the spherical coordinates r,6,cp with origin at the sphere 
centre). Energy density, e, pressures in transverse, p±, and radial, p, directions are determined 
by relations 

q(a,r) = -L- ^^T^D^ (ur), q = e,p, P± , (8.10) 

LOlX 

p(a,r) = e-2p ± , (8.11) 



where the following notations are introduced 



V i W-\ l(l + l)(2l + l)jf(y)/y*, q = p ± V' 12 * 



In the sum ( |8. 10|) u takes discrete set of values determined by the equations (8/7) and ( |8. 8| ) . 
The relation (|8.1l| ) corresponds to the zero trace of the EMT. 

The vev ( |8.10| ) are infinite. The renormalization of (0|Tjfc|0) in flat spacetime is affected by 
subtracting from this quantity its singular part (0|Tjfc|0), which is precisely the value it would 
have if the boundary were absent. Here |0) is the amplitude for Minkowski vacuum state. To 
evaluate the finite difference between these two infinities we will introduce a cutoff function 
Tp^(uj), which decreases with increasing to and satisfies the condition Tp^iuj) — ► 1, fj, — ► 0, and 
makes the sums finite. After subtracting we will allow fj, — ► and will show that the result does 
not depend on the form of cutoff: 

reg(0|T ifc |0> = lim [<0|r«|0) - <0|T tt |0)] . (8.13) 
fj,—t0 

Hence we consider the following finite quantities 

-. oo oo 1 

^' a ' r) = E4 )3T ^^)^(iS/«W 9) (iS^), x = r/a, (8.14) 

47F a i=lfc=lA=0 



where u = jH/a are solutions to the eigenvalue equations and (8^) for A = 0, 1, respec- 
tively. The summations over k in (|3D can be done by using the formula ( |3.22| ) and taking 
A = 1,B = for TM-modes (A = 0) and A = 1,B = 2 for TE-modes (A = 1) (recall that in 
( |3.18|) A^fc are zeros of J v {z) with bared quantities defined as fl3,l|) ). Note that the resulting 
sums are of type ( 3.26| ). Let us substitute in formula ( 3.22 ) 



f(z) = z 3 ifj II (z/a)D\ q, (zx), (8.15) 

with D\ q \y) defined from (jTEj) . We will assume the class of cutoff functions for which the 
function ( |8 . 1 5| ) satisfies conditions for Theorem 2, uniformly with respect to fj, (the corresponding 
restrictions for ijj^ can be easily found from these conditions using the asymptotic formulae for 
Bessel functions). Below for simplicity we will consider the functions with no poles. In this 
case ( |8.15| ) is analytic on the right-half plane of the complex variable z. The discussion on the 
conditions to cutoff functions under which the difference between divergent sum and integral 



exists and has a finite value independent any further details of cutoff function see [47]. For TE- 



and TM-modes by choosing the constants A and B as mentioned above one obtaines 

i ^ r_ r 00 , , , , , i 



iST 2 E { 2 J u 3 ^(u)Dl q) {u;r)duJ - ^ X p(z/a)l} q) (z,x)dzj , q = e,p ± ,p, (8.16) 
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where for x < 1 the functions (z, x) are defined as 

~ls 2 l+l (zx) + (l + l)sti(zx) - (21 + l)s 2 (zx)] , (8.17) 
s 2 (zx), (8.18) 



,(Px), 



,(p), 



z,x) 
z,x) 



{21 + 1 



s' t (z) 
1(1 + 1 



ZX H 

,n(P±) 



aj(z) aj(z) 



In these expressions we have introduced Ricatti-Bessel functions of imaginary argument, 

f2z 



I„(z), ei(z) 



K v (z), v = 1 + 1/2. 



1.19) 



.20) 



As (0|Tjfc|0) = lim a _ >oo (0|Tife(//, a, r) |0) the first integral in ( 8.16| ) represents the vacuum EMT 
for empty Minkowski spacetime: 



-I OO „ QQ 

47T ;=1 JO 



.21) 



This expression can be further simplified. For example in the case of the energy density one has 



r (o) 



1 

1=1 
1 

~ 2n 2 Jo 



£ / wVmH Zyf+iM + (l + l)jh x {ur) + (21 + l)jf(ur) 



du 



uJ 



^ tx (u)Y^(2l + l)jt(ur)duj 



1=0 



1.22) 



As we see the using of GAPF allows us to extract from infinite quantities the divergent part 
without specifying the form of cutoff function. Now the regularization of the EMT is equivalent 
to the omitting the first summand in ( |8.16| ), which as we saw corresponds to the contribution 
of the spacetime without boundaries. For the regularized components one obtains 



ieg q(a,r) 



1 . 00 , f'°° i \ 

' 8?r 2 Q 4 J2 J X^/a)Ff'(z,x)dz, r<a, q = e,p,pj_. 



3.23) 



By using the recurrence relations for Riccati-Bessel modified functions the expressions for the 
regularized vacuum energy density and radial pressure may be presented in the form 



■1 



8n 2 a 2 r 2 



i=i 



^2(21 + 1) dzz Xfl (z/a) 



ei(z) e[(z) 
si(z) s'^z) 



s t 2 (zx) - sf(zx) 



1 



1(1 + 1) 



(-iy z 2 x 2 



(8.24) 

where i = for q = e and i = 1 for q = p. Here we keep the cutoff factor because it plays an 
important role in the calculations of the total Casimir energy for the spherical shell (see below). 
The derivation of the vacuum densities ( |8.23 ), ( |8.24 ) given above uses GAPF to summarize 
mode sums and is based on [13, One can see that this formulae for the case of exponential 
cutoff function may be obtained also from the results of |3(| |3^], where Green function method 
is used. 

We obtained the regularized values ( 8.23| ) by introducing a cutoff function and susequent 
subtracting the contribution due to the unbounded space. The GAPF in the form ( 3.22| ) allows 
to obtain immediately this finite difference. However it should be noted that by using GAPF 
in the form ( |3.18| ) we can derive the expressions for the regularized azimuthal pressure without 
introducing any special cutoff function. To see this note that for x < 1 the function ( $.15 ) with 
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q = p± and i/j^ = 1 satisfies conditions to Theorem 2. It follows from here that we can apply 
the formula ( 3.18j ) directly to the corresponding sum over u in (8.10) or over k in (8.14) (with 
tpu = 1) without introducing the cutoff function. This immediately yields to the formula ( |8.23[ ) 
for q = p± with = 1. 

Let us consider the behaviour of the functions F^ q \z,x) in various limiting cases. By using 
the corresponding formulae for Bessel functions one obtains: 

(a) When / is fixed and z — > 



Fj e \z,x) ~ tt(1 + l/2)x 2 ('- 1 ), F[ p± \z,x) ~ nlx 2 «-V/2. 
(b) When I is fixed and zx is large 



f/ £) ~ 2F l (P± > ~ P(z + if (2i + 1; 



n(Px) 



7T 



2(zx)' 



■ exp[— 2z(l — x)], 



5.25) 



3.26) 



and the integral over z in (|3,23|) converges for all x < 1; 



(c) For large / by using the uniform asymptotic expansions for Bessel functions [21] one finds 
Fj- q) (vz,x) ~ $} 9) (z,x)exp{-2i/[77(z) - r)(zx)\] , i/ = Z + l/2 (8.27) 



with 



7Tf 



t(zx)t 3 (z) { 1 



1 
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t{z){t 2 {z) + 3) + t(zx)(5t 2 (zx) - 9) 



7T 



^t 2 (zx)t 3 (z), 



where the standard notations are used: 



t(z) 



rj(z) = yl + z 2 + In 



l + Vi-T^ 2 ' 



3.28) 
3.29) 

3.30) 



From these asymptotic formulae it follows that in 
x < 1 and the cutoff may be removed by putting Xfi = 



;T23|) and 
0. 



>.24| ) the rhs is finite for 



1,/i — > 0. From here it is obvious the 
independence of the regularized quantities on the specific form of the cutoff, on class of functions 
for which ( |3.15 ) satisfy conditions for (|3.22j ). 

From ( pT27| ), ( ^28|) and ( |8^29| ) it follows that the vev of the EMT diverge at sphere surface, 
x — > 1, due to the contribution of large I (note that, as it follows from ( B.28| ) and ( 8.2S| ), in ( 8.23 ) 
with Xa = 1 the integral over z converges at x = 1). The corresponding asymptotic behaviour 
can be found by using the uniform asymptotic expansions given above and the leading terms 
have the form 

e ~ 2 Pl _ ~ 2 rg, p ~ 2 2 , 72 • (8.31) 

307r z a(a — rj^ b07r z a z (a — rj z 

These surface divergences originate in the unphysical nature of perfect conductor boundary 
conditions and are well known in quantum field theory with boundaries. They are investigated 
in detail for various types of fields and general shape of the boundary [45, fl6f| . Eqs. ( S . 3 1| ) 



arc 



particular cases of the asymptotic expansions for EMT vev near the smooth boundary given in 
these papers. In reality the expectation values for the EMT components will attain a limiting 
value on the conductor surface, which will depend on the molecular details of the conductor. 
From the asymptotic expansions given above it follows that the main contributions to q(r) are 
due to the frequencies u> < (a — r) _1 . Hence we expect that the formulae ( 8.23 ) are valid for 
real conductors up to distances r for which (a — r) _1 -C t<Jo, with ujq being the characteristic 
frequency, such that for to > u>q the conditions for perfect conductivity are failed. 
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At the sphere centre in ( 8.23 ) I = 1 multipole contributes only and we obtain |34], |l4| 



£ (0) 



1 



2vr 2 a 4 



dzz- 



z - 1 

7+V 



2z 



+ 1 



Z + l 



2z 



Z 1 + Z + l 



-0.0381a" 



p(0) = p± (0) =e(0)/3. 



.32) 



At centre the equation of state for the electromagnetic vacuum is the same as that for blackbody 
radiation. Note that the corresponding results obtained using the uniform asymptotic expansions 



for Bessel functions [36 



] are in good agreement with 
The components of the regularized EMT satisfy continuity equation T- 
spherical geometry takes the form 



i\k 



0, which for the 



p'(r) + -(p — p±) = 0. 



.33) 



From here by using the zero trace condition the following integral relations may be obtained 



P{r) = 4 f e(t)t 2 dt = — [ r p±(t)tdt, 
r d Jo r z Jo 



.34) 



where the integration constant is determined from the relations fl8.32Q at the sphere centre. It 
follows from the first relation that the total energy within a sphere with radius r is equal to 



E(r) = 4vr / e(t)t 2 dt = 3V(r)p(r), 



.35) 



where V(r) is the corresponding volume. 

The distribution for the vacuum energy density and pressures inside the perfectly conducting 
sphere can be obtained from the results of the numerical calculations given in |?(], 37]. In their 



calculations Brevik and Kolbenstvedt use the uniform asymptotic expansions of Ricatti-Bessel 
functions for large values of order. In [13, 14] (see also [0]) the corresponding quantities are 
calculated on the base of the exact relations for these functions and the accuracy of the numerical 
results in |3(], 37 1 is estimated (~ 5%). The simple approximataion formulae are presented with 
the same accuracy as asymptotic expressions. Note that inside the sphere all quantities e, p, p± 
are negative and corresponding vacuum forces tend to contract sphere. 



9 Electromagnetic vacuum EMT outside a spherical shell 

Now let us consider the electromagnetic vacuum in the region outside of a perfectly conducting 
sphere. To deal with discrete modes we firstly consider vacuum fields in the region between two 
cocentric conducting spherical shells with radii a and b, a < b. Letting b — > oo we will obtain 
from here the result for the region under question. 

By using Coulomb gauge the complete set of solutions to the Maxwell equations can be 
written in the form similar to (Is3 



- e lwt o , u \) vgoiiua^rfXim if A = 

^ A (M) = ^(a,^) Vx[ffiiKwr)Xim] if A = 1 , (9.1) 



where as above the values A = and A = 1 correspond to the waves of magnetic (TE-modes) 
and electric (TM-modes) type, 

/ \ I ji(v)nr(x) — ii(x)ni(y) if A = ,„ „. 

Ji(y)[xni(x)\ - [xji{x)\ ni(y), if A = l 
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with ni(x) being Neumann spherical function. From the standard boundary conditions at sur- 
faces r = a and r = b one finds that possible energy levels of photon are solutions to the following 
equations 



) [ r 9\i(ua,ujr)} r=b = 0, A = 0,1. 



(9.3) 



All roots of these equations are real and simple [21]. 

The coefficients (3\i in fl9.1| ) are determined from the normalization condition ( |8.4| ), where 
now the integration goes over the region between spherical shells, a < r < b. By using the 
standard relations for spherical Bessel functions they can be presented in the form 



0X1 



UJCL 



a Jl 



ua) 



bjf(ub) 



1 \ b[ujaji(uja)} 



[a [iobji(Lob)} 



A = 
1(1 + 1) 



2 b 2 



1 + 



1(1 + 1) 

uj 2 a 2 



A = 1. 



(9.4) 
(9.5) 



From ( |3.1| ) with the electromagnetic field EMT and functions (9.1) as a complete set of solutions 
one obtains the vev in the form (^^) with 



q(a,b,r) = ^J2(2l + l)w 4 /3f z /|f (ua,ur), q = e, p, p± 



(9.6) 



UllX 



where the frequencies oj are solutions to the equations (|9.3j ), and 



f^(ua,ur) 
f^\ua,ur) 



1 + 



1(1 + 1) 



ui 2 r 2 



9\l(ua,u)r) + -g-j 



d 



d(u)r] 



(wrg\i(uja,Lur)) 



1(1 + 1) 
iv 2 r 2 



9xi(ua,u)r). 



(9.7) 
(9.8) 



It is easy to see that the eigenvalue equations (|9.3| ) can be written in terms of the function C^ B , 



defined by (41), as 



C^ B (r},ua) = 0, v = 1 + 1/2, 7] = b/a, A = 1/(1 + A), B = A, A = 0,1. 



(9.9) 



By this choice of constants A and B the normalization coefficients Q9.4| ) and ( |9.5D are related 
with the function T^ B from (POT) as: 



0X1 



(rj, ua). 



(9.10) 



This allows to use the formulae from section 4 for the summation over eigenmodes. 

As above to regularize the infinite quantities ( |9.6D we introduce a cutoff function ^^(ui) and 
consider the difference 



reg(0|T ifc |0) = lim 



(0\T ik (fi,a,b)\0) - lim lim (0\T ik (^ a, b)\0) 



(9.11) 



This procedure is equivalent to the subtraction of Minkowskian part without boundaries. 
Hence instead of fl9.6|) we consider the finite quantities 



oo 1 



E^ + ^E^^feiS^tTS/al/lftTlHJ), q = e,p x , (9.12) 



8vra 4 



(?)/„,(A) (X) r 



k=lX=0 
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where x = r/a, and toa = 7^ are solutions to the equations ( p.3[ ) or 
will use the formula ( 4.13| ) with 

h(z) = z A ij) fl {z/a)f i ${z,zx), 



To sum over k we 



(9.13) 



assuming a class of cutoff functions for which ( |9.13| ) satisfies to the conditions fl4.4|) and (|4.11|) 
uniformly with respect to jj,. The corresponding restrictions on ip„ can be obtained using the 
asymptotic formulae for Bessel functions. From ( |9.12| ) by applying to the sum over k the formula 
( |4 . 1 3| ) for the EMT components one obtains 



1 x 1 ( re 



3 / / / \ /a/ ( z > zx ) 



+- 



ef\rjz) zx^(z/a)F ( x f(z,zx) 



dz+ 



x 2 Jo 4 X \z) l(d/dy)*G xl (z,y)] 



-dz 



where we use the notations 



(») = ( 



-V 

dy) 



y=zr\ 



dy 



si(y) 



for the Riccati-Bessel functions derivatives, 

tti\i(z) = 



jf(z)+n 2 (z), 
[zUz^ + izmiz)}' 2 , 



A = 
A = 1 



and 



G xl (x,y) = e\ X \x)s l (y)-e l (y)s\ / "(x), A = 0, 1 



(A), 



F$(z,y) 



d „ , j 2 r/(/ + i) 



dy 



G\i(z,y) 



+ 



G 2 xi(z,y), 



l(l + l)G 2 xl (z,y)/y 2 . 



(9.14) 
(9.15) 

(9.16) 

(9.17) 
(9.18) 
(9.19) 



The function X/j, is determined by ( |3.18| ). To obtain the vev for the EMT components outside of 
a single conducting spherical shell with radius a let us consider the limit b —* 00. In this limit the 
second integral on the right of formula (p. 14 ) tends to zero (for large 77 = b/a the subintegrand 
is proportional to e~ 2r,z ), whereas the first one does not depend on b. Hence one obtains 



1 00 rco 
07T a 1=1 A=0,1 



f$\z,zx) 



dz, 



e,p±. 



(9.20) 



To regularize the expressions ( p.20| ) we have to subtract the Minkowskian part, namely the 
expression ( 8.21] ) . It can be easily seen that 



m\z,z X ) 



D\*\zx) 



(mg) 



m=l,2 



Here the functions D^ mq \y) are obtained from the relations ( 8.12 ) by replacing ji — ,1. 



(9.21) 



h\ m \ with 



m 



1, 2 being spherical Hankel functions, and 



j^/h^iz), A = 
[zHz)]' /[zht\z)]' , A = l 



(9.22) 
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The function (z)k {hY 1 (z)) has no zeros for < argz < tt/2 {—it/2 < argz < 0) and from 
this it follows that 



,(2), 



m=l,2 



I z z ^ li {z/a)^\z)D [ [ aq) {zx)dz = 2 / z 3 Re ^^{iz/a)^ {iz)D\ lq) {izx) dz. (9.23) 



Jmq) 



(!), 



By introducing Ricatti-Bessel functions ( |8.20| ), for the regularized components of the vacuum 
EMT outside the sphere we find 



q(a,r) 



-j^ OO POO 

-p^iE J Xn(z/a)Fl Q) (z,x)dz, 



r > a, q = e,p±,p, 



(9.24) 



where for x > 1 the functions i 7 ) (z, x) are defined as 



l=i 
(?) 



F«(*,x) = - 
F/ Px) (z,x) = (2/ + 1 



e{(z) 
Z(Z + 1 



lef +1 (zx) + + l)ef_ 1 (zx) - (21 + l)ef (zx) (9.25) 



sijz) s[{z) 
ei(z) e[{z) 



ef(zx), F l 



(p) 



2F, 



(p±) 



(9.26) 



The exterior mode sum consideration given in this section follows (T^, [Tf|]. For the case of 
exponential cutoff function the formulae fl9.24j ) and ( 9,26| ) can be obtained also from the results 
[ {3G| , 37], where Green's function formalism was used. Note that the expressions for the exterior 
components are obtained from the interior ones replacing — ► ij, ij — ► sj. In particular, the 
exterior components can be presented in the form analog to (|8.24j ). 

Let us consider the behaviour of the functions F^ q \z,x) in various limiting cases. By using 
the corresponding formulae for Bessel functions one obtains: 

(a) When I is fixed and z — > 



fI £) (z, x) tt(Z + l/2)x- 2 ^ , F^\z, x) ~ -vr(Z + 1) X - 2 ^ /2. 



(b) When I is fixed and z is large 



-I 2 (I + 1) 2 (21 + 1) 



7T 



2(zx)' 



exp[— 2z(x — 1)]; 



(9.27) 



(9.28) 



(c) For large I by using the uniform asymptotic expansions for Bessel functions [21] one finds 
Fj; q \vz,x) ~ ${ 9) (z,x)exp{-2i/[77(zx) - t?(z)]} , i/ = Z + 1/2 (9.29) 

with 



<5>f\z,x) 
^ p \z,x) 



^t{zx)t 3 (z) (l + — \t(z)(t 2 (z) +3) +t(zx)(bt 2 (zx) - 9) L (9.30) 



12zy 



71 t 2 (zx)t 3 (z), 



2x 3 



(9.31) 



whith notations ( |3.30D . 

It follows from here that for the values x > 1 the expressions (9.24) are finite and hence 
cutoff may be removed. In this case the independence of the result on specific form of cutting 
function is obvious. 

The expressions ( 9.24| ) with ip„ = 1 diverge at sphere surface. The leading terms of these 
divergences may be found using (|9.29|) and are as following 



e ~ 2p± 



1 



30ir 2 a(r — a) 3 



, p(r) 



■1 



607r 2 a 2 (r — a) 2 



(9.32) 
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and 



]im(e/p±)' = lim(p/pj 



1. 



(9.33) 



Comparing fl9~3p with flOl]) we see that the cancellation of interior and exterior leading di- 
vergent terms occurs in calculating the total energy and force acting on sphere. The same 
cancellations take place for the next subleading divergent terms as well (see below). Formulas 
( p. 32 ) are particular cases of general asymptotic expansions of the vacuum EMT components 
for conformally invariant fields near an arbitrary smooth boundary given in [45]. 
For distances far from the sphere one finds 



P± 



47r 2 a 4 x 7 



z 2 e\ (z)dz 



5a 3 



167r 2 r 7 



-4p • 



2vr 2 r 7: 



r»a. 



(9.34) 



The results of numerical calculations of the vacuum EMT components outside the sphere are 
given in [36, 37, 14]. In [36, [37]] calculations are carried out by using the uniform asymptotic 



expansions for Riccati-Bessel functions. The accuracy of this approximation is estimated in [14], 
where exact relations are used in numerical calculations. The simple approximating formulas 
with the same accuracy as those for the asymptotic calculations are presented as well. The energy 
density and azimuthal pressure are positive, and radial pressure is negative. The latter means 
that the exterior vacuum forces tend to expand sphere. As we will see below this dominates the 
interior contraction force. 

Note that the continuity equation ( |B.33 ) now may be written in the following integral form 



p(r) = -3 / e{t)t 2 dt = 4 / P±(t)tdt, 



(9.35) 



where the integration constant is determined from the asymptotic relations ( 9.32| ). From ( 8.34 ) 
and ( |9.35| ) it follows that 



E(a) = / e(r)dV = 4?ra 3 [p(a-) -p(a+)], 



(9.36) 



where E(a) is the total vacuum energy for a spherical shell with radius a, p(a±) = lim r ^o p(a±r). 
By using the expressions for p(r) given above one can obtain the following formula for the total 
energy (the same result can be obtained also by integrating the energy density) 



E(a) 



-1 OO 

27Ta t! 



dzxn{z/a)z (In \si(z)ei(z)\)' 



1 + 



1(1 + 1) 



1 



si(z)ei(z) 



1 < ^~ > P OO rJ 

Z-YpL + 1) j Q dzx^zMz-gln^-MzMz)}' 2 }. 



(9.37) 



By taking the cutting function tp^iuj) = e ^ one obtains the expression for the Casimir energy 
of the sphere derived in [33] by Green function method. Note that in this method the factor 

ip^iz/a) 



appears automatically as a result of the point splitting procedure. The 
evaluation of (9.37) leads to the result E = 0.092353/2a for the Casimir energy of a spherical 
conducting shell l30|, 31, 32, 4C, 38]. This corresponds to the repulsive vacuum force on 



the sphere. Here the cancellation of interior and exterior divergent terms in the energy density 
occurs. The discussion on cancellations of divergences between interior and exterior modes see 

M,M, MM. 
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10 Electromagnetic vacuum in spherical layer between perfectly 
conducting surfaces 

Electromagnetic vev of the EMT in the region between two cocentric perfectly conducting sur- 
faces with radii a and b, a < b, may be obtained from the results of previous section. The 
corresponding nonrenormalized components are given by ( |9.14| ), Using this formula they can be 
presented in the form 

q(a,b, r) = q(a,r) + q^'(r), a < r < b, q = e,p±,p, p = e — 2p±, (10.1) 

where q(a,r) is given by (|9.20| ), and [|l(| 



( ab )( r \ — VrtJ^n V / z^ fl (z/a)n xl (z,7 1 )F i x f(z,zx)dz, x = r/a. (10.2) 



^^E^ + i) E I 

07F a r 1=1 A=0,1 J ° 

Here the functions F^f are defined by relations (|9.18|) , (|9.19 ), n = b/a, and 



O ( \ e {X \z V )/e {X \z) 

tou(z,r,)= ( ' (10.3) 

(see notation ( |9.15| )). In ( |10.1| ) the dependence on b is contained in the summand q^ only. 
This quantity is finite for a < r < b and the regularization of q(a, b, r) is equivalent to the 
renormalization of the first summand. This procedure have been done in previous section, where 
we have seen that q(a,r) (see expressions ( |9.24D and fl9.25|) ) coincides with the corresponding 
quantity for the exterior region of a single shell with radius a. The expressions ( 10. 2| ) for 
a < r < b are finite when \i — ► and hence for these values the cuttoff function may be removed 
putting Xfi = 1- 

It can be seen that the quantities ( |10.1| ) may be written also in the form 

q(a,b,r) = q(b,r) +q (ab) (r), q = £,p ± ,p, (10.4) 

where 

9 (ab) W = ^ W E(^ + 1 ) E I zn xl (z,a)F^\z,zy)dz, (10.5) 

i— 1 A— 0,1 

with y = r/b, a = a/b, and 

"Al(*, ct ) = -73^ fl^^'n}^ m • (10-6) 

4 X \za) S ( l X) (z)-e < i X \z)s i l X) (za) 1 ' 

In ( 10. 1 ) g( ab )(r) — ► when a — > and q(b,r) coincides with the corresponding quantities inside 
a single conducting shell with radius b (the latter can be seen also by direct evaluation of q(b, r)). 
Note that in (10.5) the sum and integral are convergent for a < r <b. 

As we said above from the expressions q(a, b, r) in limiting cases a — > or b — > oo may 
be obtained the vacuum stress inside and outside of a single shell. Consider now the another 
limiting case: h = b — a = const, b — ► oo. For a/b — > 1 the main contribution in ( |10.2| ) is due 
to large I. This allows us to use asymptotic formulae for Bessel functions. For instance, in the 
case of the energy density one has 

pM*, (10.7) 
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where 

A = ^ v 3 [ e Mv(z)-v(™/b)] 



-i 



b 1 



s 3 ds 



4^4 



7T 4 6' 



16/i 4 (l + z 2 ) 2 Jo e s -l 240/i 4 (l + z 2 ) 2 ' 



(10.1 



By substituting this into (|10.7|) we receive the standard result for the Casimir parallel plate 
configuration: e = — 7r 2 /720/i 4 . 



(10.9) 

(10.10) 
(10.11) 



Let us present the quantites q = e,p,p± in the form 

q = q(a, r) + q(b, r) + Aq(a, b,r), a < r < b, 

where "interference" term may be written in two ways 

Aq(a,b,r) = qW (a,b,r) - q(b,r) 
Aq(a,b,r) = cf ah) (a,b,r) - q(a,r). 

Here q( ab ^ and q^ ab ^ are defined by relations ( |10.2| ) and ( |10.5| ). It can be seen that Aq(a,b,r) 
is finite for all a < r < b, a < b. Near the surface r = a it is convenient to use ( |10.1C| ), as for 
r — > a both summands in this formula are finite. For the same reason the formula ( |10.11 ) is 
convenient for calculations near the surface r = b. 

So far we have considered the electromagnetic vacuum in the region between two perfectly 
conducting spherical surfaces. Consider now a system consisting two cocentric thin spherical 
shells with radii a and b, a < b. In this case the vev for the EMT components may be written 
in the form 

q{a,b,r) = q(a,r)6(a - r) + q(b,r)B(r - b) + \q(a, r) + q {ab) {r) \ 9(r - a)6(b - r), (10.12) 

where 9{x) is the unit step function. By using the continuity equation ( |3.33| ) it is easy to see 
that the total Casimir energy for the system under consideration can be presentaed in the form 



= E(a) + E(b) + 4vr Itfp^ib) - a 3 p^ ab \a) 



(10.13) 



where E(i) is the Casimir energy for a single sphere with radius i, i = a,b. As it follows from 
( |i02l) and (|ia5|) the additional vacuum pressures on the spheres are equal to jig , [l2|1 



p( ab \a) 



p (ab \b) 



1 re 



dzz 



5^(2/ + 1) / dzz 
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1(1 + 1) 
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+ 1 



+ 1 



^il(z,v) ~ tt i(z,r]) 
ttu(z,cr) - h i(z,a) 



(10.14) 
(10.15) 



where and are defined by relations ( |10.3|) and ( [LP. 6 ). 

The vacuum force per unit area of the inner sphere is equal to 



(10.16) 



■! ~r i-^-L' ~ ' , ^ - -p^>(a) 

where is the force per unit area of a single sphere with radius a, and AF^ is due to the 
existence of the second sphere ("interaction" force). By similar way vacuum force acting on per 
unit area of outer sphere is 

F (b) = F (b) + Ap (b)^ Ap (b) = p{ab)( b y ( 10 17 ) 

The results of numerical calculations of quantities Aq(a,b,r), q = e,p,p±, as well as those for 
are presented in |l6|, |l2| . Note that as it follows from the results of these calculations the 
quantities ( |10. 14 ) and (10.15) are always negative, and therefore the interaction forces between 
two spheres are always attractive (as in the parallel plate configuration). The total Casimir 
energy is positive for small values of a/b and is negative for values close to 1. At a/b ~ 0.7 this 
energy is zero. 
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11 EMT vev inside a perfectly conducting cylindrical shell 



In this and next sections we will consider the case of perfectly conducting cylindrically symmetric 
boundaries. The Casimir effect for a perfectly conducting cylindrical shell was considered in 
[fiij (see also |50tl ) and for a dielectric cylinder in [51| by using the Green function formalism. 



Recently the problem is reconsidered in [52, 53 1 using the mode summation technique and in [42 



p4| , 55 1, within the framework of the zeta- function regularization scheme. In these papers global 
quantities, such as the total enehgy and stress on a shell, are investigated. Local characteristics 
of the electromagnetic vacuum are considered in 17] for the interior and exterior regions of a 



conducting cylindrical shell, and in [18| for two coaxial shells. In this papers the mode summation 
method is used combined with generalized Abel-Plana formula. Our consideration below is based 
on these works (see also [0]). The vev of the EMT for electromagnetic field inside a perfectly 
conducting cylindrical surface with radius a can be found by the way similar to the spherical 
case. As an eigenfunctions we use the vector potentials corresponding to the cylindrical waves 
of magnetic (A = 0) and electric (A = 1) type: 



P. 



-e 3 x V 4 { J m {ir) exp [i(mip + kz - ut)]} , 
(l/iu) [e 3 + (ifc/7 2 )V t ] J m {ir) exp [i(mtp + kz 



ut)] 



A 
A 



(11.1) 



where the cylindrical coordinates (r, (p, z) are used with unit vectors e^, 7 2 = uj 2 — k 2 , m is 
an integer, Vt is the transverse to the z axis part of the nabla operator. From the standard 
boundary conditions we obtain the following equations for the possible values of the quantum 
number 7: 



Jm{ia) 



0. 
0, 



A = 0, 
A = 1. 



(11.2) 



The constants (3\ m are determined from the normalization condition and are equal to 



R 2 



7 



J'mha) + (1 - m 2 /-f 2 a 2 )J'^(ja) 



TTLOa 



-r m (7a), 



(11.3) 



where T m (z) is defined by (|3.8|). As independent quantum numbers we will choose the set 
a = (mk'yX). In this case uj 2 = 7 2 + k 2 and 7 takes discrete values being solutions to ( 11. 2| ). By 
using the standard formula (|3.1|) with (11.1) one obtains 



(0|r||0) =diag(e, - Pl) -p 2 , -p 3 ), 



(11.4) 



where the energy density e, the pressure pi in direction can be presented in the form (below 
the index c will specify quantities for the cylindrical geometry) 



with 



q c (a,r) 

&\v) = 



1 

8^ 



+00 
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dk E Plmfm M , q c = e, pi , p 2 

A,7 



(11.5) 
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m 2 t2 



J m (y) + -r J m(y) 



+ Jl(y) 



m 



JM- tt + (-!)' \ J m(y) 
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(11.6) 
(11.7) 
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and p% = e—p\ —p2- The latter corresponds to the zero trace of the vacuum EMT. The quantities 
( |11,5| ) are divergent. To make them finite we introduce the cutoff function ^(7) and consider 
the finite quantities 



1c 



8vr 2 a 4 
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dk 



1 00 . W | / .(A) , s 
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k 2 a 2 + j 
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(AW ) f(l)(nW T ] 



q = e,p!,p 2 (Hi 



with x = r/a, 7a = are the roots of the equations ( pL 1 . 2|) f or A = and A = 1, corre- 
spondingly. To calculate the sums over zeros of the functions ( |11.2| ) here we use the summation 



formula obtained in section 2, namely the formula (3.21). Let us choose as a function f(z) in 
GAPF 



/(*) 



\/ z 2 + k 2 a 2 



Mz/a)f${zx) 



(11.9) 



This function has branch point on the imaginary axis and we have to use the version ( 3.21 ) with 



lower sign. Here we will assume a class of cutoff functions for which (11.9) satisfies conditions 
( |3.4p and (3.11) uniformly with respect to fi. By using the asymptotic formulae for Bessel 
functions these conditions can be easily translated in terms of tp^. We choose A = 0, B = 1 in 
the case A = and A = 1, B = in the case A = 1 (see (|3.1|)), and v = m. Using the relation 



ye 



-m/2\ 



ye 



(11.10) 



we see that the subintegrand of the first integral on rhs of the formula ( |3.21 ) is proportianal to 
ip^iz/a) — ip[j,(—iz/a). Consequently after removing the cutoff (ip^ — > 1) the contribution of the 
first integral will be zero. For this reason we shall write only the second integral on the right 
of ( |3.21| ). For the simplicity we will assume also that the cutoff function has no poles in the 
right-half plane. In this case the residue terms are zero. It can be seen that the residue term on 
the right vanishes as well. Hence by applying GAPF to the sums over zeros of Bessel functions 
in (11.8) and omitting the term which will vanish after the cutoff removing we obtain 
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(11.11) 



where the function x^iu) ls defined in ( |8.18 ). The second integral on the right of this formula 
vanishes in the limit a — > 00, whereas the first one does not depend on a. It follows from here 
that the latter corresponds to the Minkowskian part without boundaries. This can be seen also 
directly by explicit summation over m using the formula X)m=-oo ^n±m( z ) = 1- F° r instance, in 
the case of the energy density one has 
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(11.12) 



with lo 2 = z 2 + k 2 . Hence GAPF allows us to extract the contribution of unbounded space 
without specifying the cutoff function. The remained part is finite for x < 1 and fx — ► 0, and 
can be written in the form 
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2?r 3 a 4 
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dt 
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(11.13) 
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with z 2 = t 2 + y 2 and t = ka. Here we have introduced a new integration variable y and the 
prime on the summation sign indicates that the m = term is to be halved. For q = e from 
( |11.6| ) one has 



ri z 2 f$(zxe™/ 2 ) = z 2 I 2 m {zx) + (t 2 - y 2 ) 



m 



I m (zx) + 2 „I m (zx) 

Z A X Z 



(11.14) 



and it can be easily seen that the contribution of the summand containing t 2 — y 2 in Qll.13 ) is 
zero. Introducing the polar coordinates (z, 0) on the plane (t, y) for the EMT components inside 
a perfectly conducting cylindrical surface from (11.13) one finds (T] 



Qc(a,r) 



47r 2 Q 4 X/ ' J dz X^(z/a)F^>(z,x), r<a, q = e,pi, (11.15) 
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where the following notations are introduced 
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(11.17) 



In particular we see that inside the cylinder e = —P3. This relation is the same as in the case of 
the Minkowski vacuum. This is natural, as we have no constraint on z direction. On cylinder 
axis (x = 0) the m = term contributes only and we have 



e(0)=pi(0)=p2(0) 
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(11.18) 



with q'(0) = 0. The vacuum EMT satisfy continuity equation which can be written now as 

^ + -(Pi~e) = 0, (11.19) 
ar r 

or in the integral form 

E c (r) =2ir [ e(t)tdt = 7rr 2 Pl (r), r < a. (11.20) 
Jo 



To determine the integration constant here we have used the relations flll.18] ) between the EMT 
components on the cylinder axis. As we see the total energy per unit length inside the cylinder 
with radius r is equal to the radial pressure on the surface of this cylinder multiplied by the 
corresponding volume. 

Let us consider the behavior of the functions (|ll.l7j ) in two limiting cases: 

1) For fixed m and large zx from the asymptotic expansions of Bessel functions we find 



^ ! -2*(l-x) pipy) 

cm n cm n ' cm 
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2X 2 ' 
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(11.21) 



2) For large m by using the uniform asymptotic expansions of Bessel functions one obtains 
F®(mz,x) ~ ^(z,x)e X p{-2m[rj(z)-rj(zx)}}, (11.22) 



39 



with 



ffl. 
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z 5 t(zx)t 3 (z) I 1 



12m 



-z 5 t 2 (zx)t 3 (z) 



t(z)(t 2 {z) - 3) + t(zx)(5i 2 (zx) - 3) J j(11.23) 

(11.24) 



where the standard notations are used. 

It follows from here that at cylinder surface, r — > a, the expressions for the EMT components 
are divergent and near the surface the corresponding quantities are dominated by large m. 
Hence to obtain the asymptotic behaviour we can use the corresponding asymptotic formulae 
for modified Bessel functions. Then after the elementary summation over m we find the following 
asymptotic behaviour 
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(11.25) 



This formulae are special cases of the general expansions for the EMT near a smooth boundary 
of arbitrary shape [45]. Note that, as it follows from ( |11.21| ) now, unlike the spherical case, in 
11.15j ) with Xu = 1 the integral over z diverges for x = 1. 



The results of the numerical calculations for vacuum EMT components (11.15) are presented 
in [17, 12]. Note that £,pi < 0, i = 1,2 everywhere inside the cylinder. The ratio of the energy 
density to the azimuthal pressure is a decreasing function on r and 0.5 < s/p2 < 1. 



12 Vacuum EMT outside a perfectly conducting cylinder 

First we consider the vev of the electromagnetic EMT in the region between two coaxial cylin- 
drical surfaces with radii a and b, a < b. The corresponding eigenfunctions have the form ( 11.1 ) 
with replacement J m {^r) — > P\ m {^a^r), where 



P\m(x,y) 



J m (y)Y m (x) - Y m (y)J m (x), A = 1 
J m (y)Y^(x) - Y m {y)J' m {x), A = 



(12.1) 



From the boundary conditions on r = a,b one obtains that the eigennumbers 7 have to be 
solutions to the following equations 



Pim{ia,ir)\ r=b 



0, A = 1 
0, A = 



r=b 



(12.2) 
(12.3) 



These equations have infinite number of simple real solutions. Now the normalization coefficients 
(3\m are in form 
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PXm 4a 4 
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(12.4) 



where z = 7a, r/ = b/a. 

From Eq.(^]l]) it follows that the vacuum EMT has diagonal form ( |11.4j ) with components 



2 +00 p+00 
q c (a,b,r) = — ^ / dk ^^Imfxli^ V"), 9, 
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(12.5) 
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where the expressions for the functions f\m^ a ->y) are obtained from ( 11. 6| ) and (11.7) replacing 
Jm{y) P\m(ja,y). 

The eigenvalue equations ( |12.2| ) and ( 12. 3| ) can be written in terms of the function Q4.1| ) as 
(^(77,76) = 0, A = X, B = 1 - X, A = 0,1 (12.6) 



(see the notation (|3.l|) ). Note that the normalization coefficients can be expressed in terms of 

(12.7) 



the function (4.6): 
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Using these relations and introducing a cutoff function ip^ the divergent quantities ( 12.5 ) can 
be written in the form of the following finite integrosums 
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where q c = e,pi and 7a = 7m,n, n = 1, 2, . . . are the solutions to the eigenvalue equations (|12.2|) , 
( |12.3| ) or ( p.2.6 ). By choosing in the formula (4.13) 
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(as noted above this formula is valid in the case when the corresponding function has branch 
point on the imaginary axis (see also Remark to the Theorem 2)) one obtains 
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Here in accordance with ( |3.1[) and (|12.6| ) 

Jm( z ) = Jm{ z ), A = 1 

J m (z) = zJ' m {z), A = 0, 



(12.10) 
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and in similar way for other Bessel functions in ( 12.10 ). To obtain the EMT components for the 
outside region of a perfectly conducting cylindrical shell we consider the limit b —* oo. It can be 
seen that the second sum on the right of (12.10) is zero in this limit and the first one does not 
depend on b. Hence for the outside region of a single cylinder we obtain 
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To regularize we subtract from these quantities the contribution of unbounded Minkowski space- 
time which can be presented in the form (see ( |11.12j )): 
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with the function f$ defined as (|TT^ ), ([1T7|) . By using the definitions of and fm it is 
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where by definition the expression for fm q \zx) is obtained from that for fm\zx) replacing 
J m (zx) — » Hm\zx), n = 1,2, and 



J m (z)/H&\z), A = l 
J' m (z)/H$\z), A = 



Hence 
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By rotating the integration contour for z by angle 7r/2 for n = 1 and by angle —tt/2 for n = 2 
(note that the function ff^(^) (Hm\z)) has no zeros for < argz < ir/2 (— vr/2 < argz < 0)) 
and introducing Bessel modified functions for the regularized components we obtain (the reg 
sign is suppressed) 
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(the definition of the functions -F C m see below). In ( |12.18 ) the integrals are convergent for x > 1 
and /u = and hence the cutoff can be removed. In this limit the first integral is zero and for 
regularized components of EMT after transformations, similar to the interior case, one obtains 



Qc(a,r) = 2 4 E / dzx^(z/a)F^(z,x), r > a, q = e, pu 
(<?), 



where for x > 1 the functions F^fn {z, x) are defined as 
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It follows from here that e = —p3- 

The asymptotic expressions for the functions ( 12.21| ) are as follows: 
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1) For fixed m and large z: 
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2) For large m from the uniform asymptotic expansions of Bessel functions one obtains 
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As in the case of the interior components of the vacuum EMT is divergent when r — ► a with 
asymptotic behaviour 



(12.25) 
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Comparing with (11.25) we see that in calculating the total energy for the infinitely thin cylin- 
drical shell the leading divergences cancel. 

The asymptotic expressions for the vev at large distances from the cylinder axis, r>n, can 
be found from ( |12.19| ) introducing new integration variable y = zx and expanding the integrands 
over 1/x. In this limit the main contribution comes from the lowest order mode with m = 
and one obtaines 



-Xa,r) 
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Here compared to the spherical case the corresponding quantities tend to zero more slowly. 
From the continuity equation for the vacuum EMT one has the following integral relation 
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(12.27) 



where E° Mt (r) is the total energy (per unit length) outside cylinder with radius r. Combining 
this relation with ( 11.20j ) for the total vacuum energy of the cylindrical shell per unit length we 
obtain 

E c = E l c n (a) + E™\a) = vra 2 [ Pl (a-) - Pl (a+)} . (12.28) 
By taking into account the corresponding expressions for the radial pressure this yields 
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(12.29) 



In the last expression integrating by part and omitting the boundary term we obtain the Casimir 
energy in the form used in numerical calculations. The corresponding results are presented in 
] . Note that in the evaluation of the Casimir energy for a perfeclty conducting cylin- 
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drical shell by Green function method to perform the complex frequency rotation procedure an 
additional cutoff function have to be introduced (see |^|)- This is related to the abovmentioned 
divergency of the integrals over z for x = 1 . 

The results of the numerical evaluations for the energy density and pressures distributions 
(formula (|12.19| )) are presented in [0, [T^| . The energy density and azimuthal pressure in the 
exterior region are always positive, and radial pressure is negative. The ratio of the energy 
density to the azimuthal pressure is decreasing function on r, and 1/3 < e/p2 < 0.5. Note that 
this ratio is continous function for all r and monotonically decreases from 1 at the cylinder axis 
to 1/3 at infinity. 



13 Vacuum EMT between two coaxial cylindrical shells 



By using the results from previous section the vev of the electromagnetic EMT in the region 
between two coaxial conducting cylindrical surfaces may be presented in the form ( |11.4j ) with 
components 

q c (a,b,r) = q c (a,r) +q { c ab) (r), a < r < b, q c = e,pi, (13.1) 
where q c (a,r) is given by ( |l2~rl ), and 
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Here the function h(z) is denned according to (|12,9| ). As we have shown the first summand on 
the right of ( |13.1 ) presents the corresponding quantity for the vacuum outside a single perfectly 
conducting cylindrical shell with radius a. As we shall see later g( ab )(r) is finite for a < r < b at 
^ = 0, and hence regularization is necessary for q c (a, r) only. This have been done in previous 
section. We have shown that result does not depend on specific form of the cutoff function and 
can be presented in the form ( 12.19| ) and ( 12.2ip . 

In ( [13.2 ) the integral over z can be presented as a sum of two integrals along segments 
(0, \ak\) and (\ak\, oo). By using the relation ( |3.20| ) and the explicit form of h(z) it is easy to see 
that the first integral will contain the cutoff function in the form \l> „(iz / a) — ipu{— iz/a) and hence 
vanishes after the cutoff removing. For this reason below we will consider the second integral 
only. After the transformations similar to those we used to obtain ( |11.15| ), the quantities q^^ 
can be written in the form 
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Here we have introduced the notation 

Qim(z,y) = 

Qom(z,y) = 



K m (z)I m (y) - I m (z)K m (y) 

K'm( z ) L m(y) - l' m (z)K m (y). 



[13.7) 



The quantities ( [13.1 ) with ( 12.19| ) and ( |13.3| ) present the regularized vev of the EMT components 
in the region between two coaxial conducting cylindrical surfaces. Let us consider the limiting 
cases of the term ( |13.3| ). First let a/r, a/b <C 1. After replacing z — » zr\ and expanding the 
subintegrand over a/r and a/b it can be seen that 



Jab) 



(r) q c (b,r), a/r, a/b <C 1, r < b, 



'13. 



where q c {b, r) are the components for the vacuum EMT inside a single cylindric shell with radius 
b (see ( |iTT5| ), ( gngP ). 

When a — > b the sum over m in ( |13.3| ) diverges. Consequently for b — a <C b the main 
contribution to is due to large m. By using the uniform asymptotic expansions for Bessel 
functions in this limit one obtaines 
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(13.9) 



which coincides with the corresponding quantity for the Casimir parallel plate configuration. 



From ( 13.1 ) and ( |13.3| ) it can be seen that the vev of EMT components can be written also 
in the form 

q c (a,b,r) = q c {b,r) +q ( f b) {r), a < r < b, q c = e, Vi . (13.10) 
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Here q c (b,r) are vev inside a single cylindrical surface with radius b (see ( 11.15 ), ( 11.16 ) with 
replacement a — ► b), and 



(r) 
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where y = r/b, a = a/b, and 
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(13.12) 
(13.13) 



The quantites ( 13.11| )are finite for all a < r < b and diverge on surface r = a. 

From the above it follows that if we present the vacuum EMT components between cylindrical 
surfaces in the form 

q c = q c (a, r) + q c (b, r) + Aq c (a, b, r) (13.14) 

then the quantities 



Aq c (a,b,r)=qi ab Hr)-q c (b,r) 



~(ab) 



(r) - q c (a,r) 



(13.15) 



are finite for all r from a < r < b. In ( 13. 15|) the first presentation is convenient near the surface 
r = a, as in this case both summands are finite. Similarly the second presentation is convenient 
near r = b. 

Let us consider a system of two coaxial thin cylindrical shells with radii a and b, a < b. The 
vacuum EMT components may be written in the form 



q c (a,r)0(a-r) + q c (b,r)0(r-b)+ q c (a,r) + qf(r) 9{r - a)0(b - r). 



13.16) 



Similar to the spherical case using the continuity equation ( |11.19| ) the total Casimir energy for 
this system may be written as 



= E c {a) + E c {b) + nb 2 p^\b) - vra 2 #' 



13.17) 



where E c (i) is the Casimir energy for a single cylindrical shell with radius i, i = a,b. For the 
additional vacuum pressures on the cylindrical surfaces from ( 13.3; ) and ( |1 3 . 1 1 ) one has: 
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13.18) 
13.19) 



where n c Xm (y,z) and n c lm (a,z) are defined in ( pl|) , (|l3j), (p!2] ) and (p!3|) . In ( pl^ ) and 
( |13.19| ) the first summands in braces come from the magnetic waves contribution, and second 
ones from the electric type waves. 

Let us now consider the interaction forces between cylindrical surfaces. The force acting per 
unit area of the inner surface can be presented in the form 



= + AF^ , AF^ = -p^ (a) , 



13.20) 



where is the force acting on a single cylindrical surface with radius a, and AF^ is additional 
force due to the existence of the outer surface and is determined from ( 13.18 ). The latter is 
finite without additional subtractions. 
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By similar way the force acting per unit area of the outer cylinder 



= F® + AFP , AF^ = (b) (13.21) 



where additional term AF^ is due to the existence of the inner cylinder and is defined by 



DTIaD - 

The results of the numerical calculations for quantities Aq(a, b, r) are given in [18|. Note that 
the sign of Ae and Ap c \ is the same as in the case of interior of the parallel plate configuration. 
In particular the additional forces AFc always have attractive nature. 



14 Summary 

In the present paper we considered a possible way for generalization of Abel-Plana summation 



formula, proposed in [11|. The generalized version contains two meromorphic functions f(z) and 
g(z) and is formulated in the form of Theorem 1. The special choice g(z) = — if(z) cot ttz with 
f(z) being an analytic function in the right half-plane gives APF with additional residue terms 
coming from the poles of f(z). Another consequence from GAPF is the summation formula 
( |2.17 ) over the points with integer values of an analytic function. An application of this formula 



to the Casimir effect is given in flOfl . 

Further we consider the applications to the series and integrals involving Bessel functions. 
First of all, in section 3 choosing the function g(z) in the form (|3.2j ) we derive two types of 
summation formulae for the series J2kTv(^v,k)f(^v,k) (the definition T u (z) see ( |3.8| )) with A^ 
being the zeros of the function J u (z) = AJ u (z) + BzJ' u {z). Such a type of series arises in a 
number of problems of mathematical physics with spherical and cylindrical symmetry. As a 
special case they include Fourier-Bessel and Dini series (see pQl ). Using the formula ( [3.18 ) the 



difference between the sum over zeros of J v {z) and corresponding integral can be presented in 
terms of an integral involving Bessel modified functions plus residue terms. For a large class of 
functions the last integral converges exponentially fast and is useful for numerical calculations. 
The mode summation method for calculating the vev of the EMT inside perfectly conducting 



spherical and cylindrical shells used in [13, 14, 17] is based on this formula. In this method 



the independence of the regularized EMT components on specific form of the cutoff function 
becomes obvious. APF is a special case of ( |3.18 ) with v = 1/2, A = 1, B = and an analytic 



function f(z). Choosing v = 1/2, A = 1, B = 2 we obtain APF in the form ( 2.15| ) useful for 



fermionic field calculations. Note that the formula ( |3.18 ) may be used also for some functions 



having poles and branch points on the imagianary axis. 

The second type of summation formulae, formula ( |3.33D , considered in subsection 3.2 (The- 
orem 3), is valid for functions satisfying condition ( |3.27| ) and presents the difference between 
the sum over zeros of J u (z) and corresponding integral in terms of residues over poles for f{z) 
in the right half-plane (including purely imaginary ones). It may be used to summarize a large 
class of series of this type in finite terms. In particular, the examples we found in literature, 
when the corresponding sum may be presented in closed form, are special cases of this formula. 
A number of new series summable by this formula and some classes of functions to which it can 
be applied is presented. 

In Section 4 we consider applications to the series of type Y^ k ^^(^ 7u,k)h(ju,k) (with 



Tjy (A, z) defined as (4.6)), where 7^ are zeros of the function J v {z)Y v (Xz) — J u {\z)Y u {z). 



The corresponding results are formulated in the form of Corollary 2 and Corollary 3. Using 



the formula ( 4.13 ) the difference between the sum and corresponding integral can be expressed 



as an integral containing Bessel modified functions plus residue terms. For the large class of 



functions h(z) this integral converges exponentially fast. The formula of the second type, (|4.16|) , 



allows to find in closed form the sums of some types of the series over 7^. To evaluate the 
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corresponding integral the formula can be used derived in section 7. This yields to the another 
summation formula, ( 4.19| ), containing residue terms only. The similar formulae can be obtained 
for the series over zeros of the function J' v {z)Y u (\z) — J v {Xz)Yl{z) as well. Note that the several 
examples we found in literature when the corresponding sum was evaluated in closed form are 
special cases of the formulae considered here. We present new examples and some classes of 
functions satisfying the corresponding conditions. The possibilities are endless. 

The results from GAPF for the integrals of type p. v. F(x)J v {x)dx (see notation (|3,l|)) 
and p.v. / °° F(x)[J v {x) cos 5 + Y v (x) sin 5]dx are considered in section 5. The corresponding 
formulae have the form (|5.4|) , ( |5.7|) and ( |5.18| ). In particaular the formula (5.4) is useful to express 
the integrals containing Bessel functions with oscillating subintegrand through the integrals of 
modified Bessel functions with exponentially fast convergence. The results obtained in [24 are 
special cases of these formulae. The illustrating examples of applications of the formulae for 
integrals are given in section 6 (see (|6.3p-(p\7|) and ( |6.13| )- ( J6.17 )). Looking the standard books 
(see, e.g., @, [^]-[^9|) one will find many particular cases which follow from these formulae. 
Many new integrals can be evaluated as well. We consider also two examples of functions having 
purely imaginary poles, ( |6.19 ) and (|6,23| ), with corresponding formulae ( |6,20 ) and ( |6.24j ) (two 
special cases of these formulae see |2C|). 



By choice of the functions f(z) and g(z) in accord with (|7.l[) formulae (|7.7|) and (|7.16|) for 
integrals of type 

J u {x)Y^\x) - J^Xx)Y u {x) 



P - V -io J3W + Y?(x) FiX)dX 

can be derived from GAPF. The corresponding results are formulated in the form of Theorem 5 
and Theorem 6 in section 7. The several examples for the integrals of this type we have been able 
to find in literature are particular cases of the formula ( [7.7[ ) . New examples when the integral 
is evaluated in finite terms are presented. Some classes of functions are distinguished to which 
the corresponding formulae may be applied. 

In the following sections, based on [|13]]-|18|], the physical applications of the summation for- 
mulae obtained from GAPF are reviewed. We consider the vacuum expectation values of the 
energy-momentum tensor for the electromagnetic field inside and outside the perfectly conduct- 
ing spherical and cylindrical shells, as well as between two conducting cocentric spherical and 
coaxial cylindrical surfaces. The corresponding mode sums contain the series over zeros of Bessel 
functions and their combinations. The application of the summation formulae from sections 3 
and 4 allows (i) to extract from corresponding divergent quantities the contribution of the un- 
bounded space in explicitly cutoff independent way, and (ii) to obtain for the regularized values 
strongly convergent integrals. To compare note that in the Green function method after the 
subtraction of the Minkowskian part the additional complex frequency rotation is used. In the 
regularization scheme based on the summation formulae of APF type the complex frequency ro- 
tation is made automatically. The corresponding global quantities such as total Casimir energy 
or forces acting on the surfaces can be obtained from the EMT components. It is shown that in 
the geometries with two surfaces the additional vacuum forces due to the existence of the second 
surface always have attractive nature. In the limiting case of the large radii the corresponding 
results for the Casimir parallel plate configuration are obtained. 

Of course the applications of the summation formulae obtained from GAPF are not restricted 
by the Casimir effect only. Similar types of series will arise in considerations of various physical 
phenomenon near the boundaries with spherical and cylindrical symmetries, for example in 
calculations of the electron self-energy and the electron anomalous magnetic moment (for the 
similar problems in the plane boundary case see, e.g., 56, 57] and references therein). The 
dependence of these quantities on boundaries originates from the modification of the photon 
propagator due to the boundary conditions imposed by the walls of the cavity. 
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